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-p  Electrode  effects  like  spot  formation  and  spot  motion  can  cause 
high  erosion  on  the  cathode  and  backplate  of  an  MPD  thruster. 
These  effects  are  qualitatively  and  partially  quantitatively 
explained  by  means  of  a  unique  theoretical  approach.  In  conjunc¬ 
tion  with  the  theoretical  work,  experiments  on  a  special  elec¬ 
trode  test  apparatus  are  conducted  in  order  to  measure  cathode 
erosion  and  to  determine  spot  formation  and  motion  and  the  at¬ 
tachment  mode  (spotty;  diffuse  thermionic)  under  various  opera¬ 
ting  conditions. 

Several  numerical  model  calculations  were  developed  which  allow 
to  determine  the  flow  arc  discharge  region  and  the  performance 
(thrust  and  specific  impulse)  of  a  nozzle  type  and  a  cylindrical 
MPD  thruster  as  a  function  of  mass  flow  rate  and  electric  cur¬ 
rent.  For  the  nozzle  type  thruster  the  heat  flux  at  the  nozzle 
throat  was  calculated  for  various  currents  and  compared  with  ex¬ 
perimental  measurements. 

'By  means  of  these  performance  calculations  the  onset  conditions 
were  calculated  based  on  the  fact  that  due  to  magnetic  contrac¬ 
tion  (pinch  effect)  the  plasma  density  becomes  zero  at  the  anode 
(onset  theory  of  Hiigel).  The  results  of  these  calculations  are 
in  excellent  agreement  with  those  of  the  experiments. 
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1 .  Summary 

The  objectives  of  this  work  are  to  investigate  the  still  unknown 
effects  on  electrodes  and  to  determine  the  acceleration  mecha¬ 
nism  and  the  stability  behavior  within  the  flow  arc  discharge 
region  of  a  cylindrical  and  a  nozzle-type  MPD  accelerator. 

Electrode  effects  like  spot  formation  and  spot  motion  can  cause 
high  erosion  on  the  cathode  and  backplate  of  an  MPD  thruster. 
These  effects  are  qualitatively  and  partially  quantitatively 
explained  by  means  of  a  unique  theoretical  approach.  In  conjunc¬ 
tion  with  the  theoretical  work,  experiments  on  a  special  elec¬ 
trode  test  apparatus  are  conducted  in  order  to  measure  cathode 
erosion  and  to  determine  spot  formation  and  motion  and  the  at¬ 
tachment  mode  (spotty;  diffuse  thermionic)  under  various  opera¬ 
ting  conditions. 

Several  numerical  model  calculations  were  developed  which  allow 
to  determine  the  flow  arc  discharge  region  and  the  performance 
(thrust  and  specific  impulse)  of  a  nozzle  type  and  a  cylindrical 
MPD  thruster  as  a  function  of  mass  flow  rate  and  electric  cur¬ 
rent.  For  the  nozzle  type  thruster  the  heat  flux  at  the  nozzle 
throat  was  calculated  for  various  currents  and  compared  with  ex¬ 
perimental  measurements. 

By  means  of  these  performance  calculations  the  onset  conditions 
were  calculated  based  on  the  fact  that  due  to  magnetic  contrac¬ 
tion  (pinch  effect)  the  plasma  density  becomes  zero  at  the  anode 
(onset  theory  of  Hiigel).  The  results  of  these  calculations  are 
in  excellent  agreement  with  those  of  the  experiments. 
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2.  Electrode  Investigation 
2. 1  General  Remarks 

The  cathode  attachment  of  an  arc  is  either  spotty  and  governed 
by  field  emission  or  (seemingly)  diffuse  and  explainable  by 
thermionic  emission.  While  the  first,  spotty  mechanism  occurs  on 
cold  cathode  surfaces  of  pulsed  MPD  thruster  devices  and  during 
the  starting  phase  of  continuously  running  devices,  the  thermi¬ 
onic  mode  prevails  when  the  cathode  surface  is  extremely  hot, 
i.e.  reaches  temperatures  above  at  least  2000°  K.  In  order  to 
sustain  such  temperatures  the  cathode  material  must  be  of  a 
highly  refractory  metal  like  tungsten  or,  what  works  best  so 
far,  an  alloy  of  thoriated  tungsten  (2%  ThC>2).  Since  the  thorium 
oxide  melts  at  about  200°  K  below  the  tungsten  melting  point  and 
any  stable,  spotty  discharge  requires  a  molten,  even  evaporating 
area,  one  may  conclude  that  also  at  so-called  thermionically 
emitting  cathodes  many  microspots  are  present  which  are  highly 
nonstationary  and  evenly  distributed  over  a  larger  hot  surface. 
This  conception  is  also  in  agreement  with  measured  erosion  rates 
of  both  attachment  modes,  both  of  which  may  be  explained  by  one 
and  the  same  spot  mechanism:  1  a)  the  diffuse  attachment  area 
which  consists  of  many  evenly  distributed  and  highly  nonstation¬ 
ary  microspots,  each  of  which  carries  a  low  current  of  a  few 
tens  of  amperes  and  b)  a  typical  spotty  attachment  presented  by 
one  or  several  hot  spots,  each  of  which  has  currents  up  to  the 
tens  and  even  hundreds  of  amperes.  These  latter  bigger  spots  are 
likely  the  result  of  clustered  microspots  which  are  somehow  hin¬ 
dered  in  spreading  over  a  larger  cathode  area.  One  goal  of  this 
investigation  was  to  understand  and  to  predict  the  motion  and 
behavior  of  these  spots.  It  was  shown1  that  the  overall  erosion 
rate  of  a  cathode  does  depend  on  the  spot  current  and  not  neces¬ 
sarily  on  the  overall  total  current.  It  would  therefore  be  ad- 
vantageous  to  have  many  microspots,  each  carrying  only  a  small 
current  and  all  spread  over  a  large  area,  than  to  have  only  one 
or  a  few  larger  ones  (or  many  clustered  smaller  ones  which  stick 
together  and  build  one  or  a  few  large  spot  sites)  of  fairly 
large  currents. 

1  H.  O.  Schrade,  M.  Auweter-Kurtz  and  H.  L.  Kurtz,  "Cathode 
Erosion  Studies  on  MPD  Thrusters",  18th  International  Electric 
Propulsion  Conference,  AIAA-85-2019,  Alexandria,  VA,  1985 
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2. 2  Theory 

Any  active  arc  spot  may  be  characterized  by  a  crater-like  pit 
with  an  inner  molten  layer,  above  which  is  located  a  high  pres¬ 
sure  plasma  cloud.  This  high  pressure  plasma  cloud,  caused  by 
ohmic  heating  and  heavy  evaporation  of  the  electron  emitting  in¬ 
ner  surface  layer,  expands  out  into  the  ambient  lower  pressure 
gas  sphere  and  forms  a  plasma  jet.  This  jet  contains  a  current 
carrying  channel  which  emanates  from  the  high  pressure  plasma 
cloud  of  the  crater  and  which  electrically  connects  the  low 
pressure  interelectrode  plasma  and  hence  the  anode  with  the 
cathode  (see  Fig.  1). 


By  means  of  a  unique  theory,  the  stability  behavior  of  such  a 
current  carrying  channel  within  an  axial  jet  has  been  discussed 
and  presented  at  the  19th  International  Electric  Propulsion  Con¬ 
ference  in  Colorado  Springs,  Colorado,  May  11-13,  1987  (see  Ap¬ 
pendix  II,  Cathode  Phenomena  in  Plasma  Thrusters).  It  turns  out 
that  under  certain  conditions  this  channel  becomes  unstable  with 
respect  to  small  deviations  off  its  straight  configuration.  The 
current  carrying  plasma  channel  then  starts  to  bend  more  and 
more  and  eventually  touches  the  cathode  surface  at  a  nearby  area 
of  the  spot  (see  Fig.  2).  This  cathode  surface  area  is  heated 
up,  melts  and  eventually  a  new  spot  is  created. 


Three  cases  shall  be  discussed:  the  spot  discharge  I)  without 
an  external  magnetic  field,  II)  with  a  magnetic  field  parallel 
to  the  cathode  surface  and  III)  in  which  the  magnetic  field  is 
inclined  with  respect  to  the  surface. 


I)  If  no  external  magnetic  field  Bu  is  present,  instability  oc¬ 
curs  when  the  dimensionless  quantity 
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Main 

Current  Carrying  Plasma  Channel 


Fig  1  Scheme  of  a  Spot  Discharge  with  a 
Straight  Discharge  Channel 


6 


decreases  below  zero.  Herein  Is  is  the  current  carried  by  the 

spot,  /  P  v  2  dA 

A  z  z 

z 

is  the  impulse  transport  along  the  channel  axis  and  f  is  a  fac¬ 
tor  of  about  0.8  which  weakly  depends  on  the  current  density 
distribution  over  a  cross  sectional  area  (Az)  of  the  discharge 
channel,  e  can  also  be  approximated  by  the  expression  (1b)  which 
depends  on  a)  the  pressure  difference  (pc  -  pm)  between  the 
high  pressure  (pc)  spot  plasma  and  the  ambient  pressure  (p*,) 
(see  Fig.  1)  and  which  depends  on  b)  the  average  magnetic  pinch 
pressure 

P0 

8  n  1  s^ 

close  to  the  crater  orifice.  The  factor  3.3  is  taken  strictly  as 
a  slowly  varying  function  depending  on  the  crater  shape  and  the 
current  density  distribution  across  the  channel  cross  section 
just  above  the  crater  orifice;  it  can  be  well  approximated  by 
3.3  (±  15%).  (See  Appendix  I.) 

In  order  to  keep  the  spot  current  down  and  hence  the  erosion 
rate,  one  therefore  should  try  to  get  a  spot  unstable  before  it 
reaches  a  high  current  rate.  According  to  our  results  this  oc¬ 
curs  earlier  when  the  ambient  pressure  pm  (see  eq.  (1b)  for  e) 
is  larger  than  when  it  is  smaller,  which  agrees  with  certain  ex¬ 
perimental  findings  (see  below). 

II)  If  a  magnetic  field  parallel  to  the  cathode  surface  is  ap¬ 
plied,  one  must  distinguish  between  three  discharge  channel  re¬ 
actions  (see  Fig.  3);  a)  With  respect  to  a  small  disturbance  S 
the  channel  adjusts  itself  to  one  with  a  straight  discharge 
axis;  this  reaction  is  called  "straight  stable",  b)  With  respect 
to  a  small  disturbance  S  the  channel  deviates  even  more  from  its 
original  straight  configuration  but  adjusts  itself  to  a  "curved- 
stable"  configuration  (i.e.  any  small  deviation  off  this  curved 
state  would  cause  the  channel  to  adjust  to  this  curved  configu¬ 
ration)  .  c)  With  respect  to  a  small  disturbance  S  the  discharge 
channel  bends  more  and  more  and  eventually  touches  the  cathode 
surface  at  a  nearby  area  as  shown  before;  this  last  channel  re¬ 
action  is  called  "unstable". 


Fit]  t  Stable  and  unstable  channel  behavior  depending  on  i  and  the 
direction  of  the  disturbance  & 
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Depending  on  the  amouunt  of  e  given  in  eqs.  (la)  and  (1b)  and  on 
the  direction  of  the  disturbance  t,  the  discharge  channel  reac¬ 
tion  can  now  vary  between  a)  a  straight  stable,  b)  a  curved 
stable  and  c)  an  unstable  one  (see  Fig.  4). 


For 


e  >  e 

cr  it 


8* 

HO  Is 


(2) 


where  a  depends  on  the  current  density  distribution  (see  Appen¬ 
dix  II)  and  is  expected  to  be  smaller  than  about  +  0.2,  the  dis¬ 
charge  channel  reaction  is  straight  stable. 


If  e  falls  in  the  range  0  <  e  <  +£crit  the  discharge  channel 
can  be  stable,  curved  stable  or  unstable  depending  on  the  direc¬ 
tion  of  the  disturbance  vector  S.  For  any  small  disturbance  in 
the  retrograde  direction,  i.e.  opposite  to  (]  x  5),  the  current 
channel  is  unstable;  for  any  small  disturbance  in  the  amperian 
direction  the  channel  will  assume  a  curved  stable  configuration. 


and  for  disturbances  parallel  or  anti-parallel  to  B0  the  channel 
reacts  in  a  stable  fashion.  The  range  of  the  azimuth  angle 
29crit  (~ecrit  <  0o  <  +ecrit'  where  0o  determines  the  di¬ 
rection  of  the  disturbance)  at  which  the  channel  still  reacts  in 
a  stable  manner  thereby  decreases  the  more,  the  smaller  e  be¬ 
comes.  Simultaneously  the  range  in  which  any  small  disturbance 
leads  to  an  unstable  or  curved  stable  channel  behavior  increases 
with  decreasing  e.  If  e  equals  zero,  any  small  disturbance  to 
the  retrograde  side  leads  to  an  unstable  behavior  of  the  dis¬ 
charge  channel,  and  any  disturbance  to  the  amperian  side  leads 
to  a  curved  stable  reaction.  If  e  decreases  further,  even  below 
zero,  the  angular  range  in  which  a  disturbance  leads  to  a  curved 
stable  channel  behavior  becomes  smaller  and  smaller  and  finally 
ceasAs  when  e  reaches  the  negative  critical  value  £Crit*  Fcr 
z  <  -Gcrit  a  small  disturbance  in  any  arbitrary  direction 
causes  an  unstable  channel  reaction. 


The  consequence  of  this  instability  behavior  of  the  discharge 
channel  is  the  fact  that  when  e  decreases  below  ecrit  the 
spots  should  preferably  jump  in  the  retrograde  direction  (see 
also  Appendix  II). 
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III)  If  the  external  magnetic  field  is  inclined  with  respect  to 
the  cathode  surface,  this  "retrograde  spot  motion"  (discontinu¬ 
ous  jumps)  occurs  preferably  under  a  certain  angle  with  respect 
to  the  retrograde  direction  (see  Appendix  II).  This  experimen- 
tally  observed  phenomenon,  known  as  Robson  drift  motion,  like 
the  phenomenon  of  retrograde  motion  is  therefore  explainable  by 
one  unique  theoretical  approach. 

In  a  future  step  the  behavior  of  the  spot  discharge  will  be  con¬ 
sidered  in  a  magnetic  field  perpendicular  to  the  surface  of  the 
cathode . 


2  A.  E.  Robson,  "The  Motion  of  an  Arc  in  a  Magnetic  Field", 
Proceedings  of  the  Fourth  International  Conference  on  Ionization 
Phenomena  in  Gases,  Vol.  lib,  p.  346,  Aug.  1959 


2.3  Experiment 
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The  goal  of  these  electrode  experiments  is  to  investigate  the 
cathode  attachment  mechanism  and  to  determine  erosion  rates 
under  various  operating  conditions  like  cathode  material,  cath¬ 
ode  temperature,  applied  magnetic  field,  ambient  pressure,  over¬ 
all  current  etc.,  and  finally  these  experiments  shall  guide  and 
prove  (or  deny)  the  validity  of  the  theory.  The  experimental 
test  rig  as  shown  in  Fig.  5  which  has  been  explained  in  the 
Final  Scientific  Report3  served  as  an  overall  investigation  on 
cold  and  hot  cathodes  of  different  materials.  The  results  so  far 
of  erosion  rate  measurements  for  different  materials  and  for 
different  current  pulse  numbers  on  cold  cathodes  are  plotted  in 
Table  1.  There  is  an  indication  that  on  polished  surfaces  the 
erosion  rates  are  slightly  smaller  than  on  unpolished  ones. 
After  several  shots  of  current  pulses  (duration  2  ms  each) , 
however,  any  polished  surface  acquires  a  "sand  blasted"  look 
which  stems  from  a  bombardment  of  arc  spots  leaving  a  field  of 
superimposed  craters  on  the  surface. 

Therefore,  as  can  be  seen  from  Table  1,  it  is  obvious  that  the 
erosion  rates  measured  between  polished  and  unpolished  materials 
differ  less,  the  higher  the  number  of  current  pulses  with  the 
same  electrode  sample.  Interesting  to  note  here  is  also  the 
somewhat  higher  erosion  rate  with  pure  tungsten  than  with  thori- 
ated  tungsten  electrodes. 

Photomicrographs  of  the  cathode  surfaces  which  were  exposed  to 
one  single  shot  showed  typical  distinct  craters  and  grooves,  the 
latter  of  which  may  be  interpreted  as  an  overlapping  string  of 
craters  (see  Figs.  6  and  7).  The  main  difference  between  both 
figures  is  the  width  or  size  of  the  grooves  and  craters.  At 
higher  ambient  pressure  (p»  =  0.5  mbar;  see  Fig.  7)  the  aver¬ 
age  size  of  the  grooves  (or  craters)  is  smaller  than  at  lower 
ambient  pressure  (p»  =  5  •  1 0~ 4  mbar;  see  Fig.  6).  One  also  ob¬ 
serves  during  the  experiments  that  at  higher  pressure  the  arc 

3  H.O.  Schrade,  M.  Auweter-Kurtz  and  H.L.  Kurtz,  "Basic  Proc¬ 
esses  of  Plasma  Propulsion",  Final  Scientific  Report,  AFOSR 
Grant  82-0298,  Jan.  1987 
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Table  1  Erosion  rates  for  different  cathode  surfaces. 


Number  of 
Current  Pulses 

Htatu re  4 

Erosion  Rate 

Thoriated  Tungsten 

300 

35  ug/C 

(2%  Th02 ) 

(not  polished) 

5000 

41  ug/C 

Thoriated  Tungsten 

100 

29  ug/C 

(2%  Th02) 

500 

34  ug/C 

( pol  ished ) 

6000 

42  ug/C 

Pure  Tungsten 

100 

44  ug/C 

(not  polished) 

500 

50  tic/ C 

5000 

54  ug/C 

Copper 

200 

73  ug/C 

(not  polished) 

800 

77  ug/C 

Copper 

100 

74  ug/c 

(polished) 

800 

75  ug/c 

Fig . 


Photomicrographs  of  the  arc  traces  left  on  the  surface 
of  a  cold  thoriated  tungsten  cathode  (rectangular  cur¬ 
rent  pulse  I  =  130  A,  electric  load  Q  =  0.26  Coul , 

p®  =  5 • 1 0~  4  mbar) . 


Fig.  7  Photomicrographs  of  the  arc  traces  left  on  the  surface 
of  a  cold  thoriated  tungsten  cathode  (rectangular  cur¬ 
rent  pulse  I  =  144  A,  electric  load  Q  =  0.29  Coul , 

p®  *  0,5  mbar )  . 
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Photomicrographs  of  the  arc  traces  left  on  the  surface 
of  a  hot  thoriated  tungsten  cathode  (T  =  2000°K,  rectan¬ 
gular  current  pulse  I  =  360  A,  electric  load  Q  =  0.72 
As,  p0  =  5 *10“ 4  mbar). 


4 


Fig .  8 


1  7 


spots  spread  more  quickly  from  the  original  impact  area  close  to 
the  anode  along  the  cathode  rod  than  it  is  observed  under  the 
low  pressure  condition  when  the  arc  sticks  more  to  the  original 
impact  area.  This  spot  behavior  is  an  indication  that  proves  the 
theoretical  result  as  discussed  before,  that  at  higher  ambient 
pressure  (p®)  e  should  become  more  quickly  negative  and  there¬ 
fore  the  discharge  channel  more  quickly  unstable  and  the  spots 
more  moveable.  Nevertheless,  this  indication  has  to  be  checked 
more  carefully  and  cannot  simply  be  extrapolated  to  pressures 
above  10  and  more  mbar. 

Experiments  with  hot  cathodes  again  show  spotty  arc  traces  (see 
Fig.  8)  but  the  craters  look  smoother  than  in  the  case  of  cold 
cathodes . 

Reproducible  and  reliable  erosion  rate  measurements  on  hot  cath¬ 
odes  failed  so  far  with  the  old  test  apparatus,  despite  several 
days  of  outgasing  of  the  tank  and  improvements  in  the  insulation 
material.  Therefore  a  new  electrode  test  apparatus  is  being  de¬ 
signed  and  built. 
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3 .  MPD  Thruster  Performance  Calculations 

In  order  to  predict  MPD  arc  thruster  performance,  several  model 
calculations  have  been  developed  and  are  explained  in  Appendix 
III,  "Numerical  Modeling  of  the  Plow  Discharge  in  MPD  Thrust¬ 
ers".4  This  paper  was  presented  at  the  19th  International  Elec¬ 
tric  Propulsion  Conference  in  May  1987.  The  following  section 
therefore  contains  only  more  recent  theoretical  results  concer¬ 
ning  the  two  MPD  thruster  types  shown  in  Figs.  9  and  10.  Both 
thrusters  are  also  experimentally  tested  in  our  institute,  al¬ 
lowing  certain  comparison  checks  between  theory  and  experiment. 

3 . 1  Nozzle  Type  Thruster 

The  flow  discharge  model  calculations  are  based  on  a  combination 
of  a  steady  state,  one-d imens ional  MPD  flow  and  a  two-dimension¬ 
al  current  coutour  line  treatment.  The  plasma  within  the  nozzle 
was  assumed  fully,  singly  ionized,  whereby  the  expansion  process 
of  the  electrons  was  taken  isothermally  and  that  of  the  ions 
adiabatically.  Heat  losses  and  friction  were  neglected.  The  cal¬ 
culations  were  carried  out  for  three  different  jet  plume  condi¬ 
tions4  as  shown  in  Figs.  11a,b,c  and  12a, b,c  and  for  the  follow¬ 
ing  two  electrode  conditions. 

One  time  the  electric  field  vector  E  and  the  other  time  the  cur¬ 
rent  density  vector  3  were  taken  perpendicularly  to  the  elec¬ 
trode  surfaces.  The  results  for  both  electrode  conditions,  how¬ 
ever,  differ  only  slightly  as  shown  in  Fig.  13a  and  b.  Both  cur¬ 
rent  contour  line  figures  are  based  on  the  hyperbolic  plume  of 
Fig .  11b. 

The  computer  code  for  these  semi-two-dimensional,  steady  state, 
nozzle  type  thruster  calculations  were  rewritten  for  the  Cray  II 
computer  which  became  available  at  the  University  of  Stuttgart 
about  a  year  ago. 


4  H.  0.  Schrade,  M.  Auweter-Kurtz  and  H.  L.  Kurtz,  "Basic 
Processes  of  Plasma  Propulsion",  Final  Scinetific  Report,  AFOSR 
Grant  82-0298,  Jan.  1987 
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Fig.  11  Different  plasma  jet  plume  conditions 
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fi 9  .  12  current  contour  lines  for  different 

expansion  models  (  I=2000A,  m*0.8g/s  ) 


a)  current  density  perpendicular  to  electrode  surface 
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13  Current  contour  lines  within  a  nozzle  type  thruster  for 
two  different  electrode  conditions. 
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Tne  nozzle  throat  conditions  are  based  on  the  premise  that  the 
arc  column  fills  the  cross-sectional  throat  area  and  that  only 
axial  components  of  the  velocity  and  current  density  vector  are 
present  ("fully  developed  arc";  see  loc.  cit.).  Neglecting  radi¬ 
ation  and  friction,  the  energy  equation  reduces  then  to 


1 

r  dr 


(r  X 


dT  -i 
d?J 


+ 


1 

r  dr 


(r  qr) 


(3) 


i.e.  within  the  nozzle  throat  the  radial  heat  conduction  loss  is 
balanced  by  ohmic  heating.  qr  is  the  radial  heat  flux  vector, 
X  the  heat  conduction  coefficient  and  a  the  electric  conducti¬ 
vity.  The  first  integration  of  eq.  (3)  results  in  the  heat  flux 
per  unit  axis  length  given  by 

r o  ,2 

2 nr o  q„(ro)  =  2n  J  -j—  r  dr  (4) 


By  modeling  the  current  density  distribution  across  the  throat 
area  (nro2)  through  a  paraboloid  of  grade  n  as 

j  =  j  (i  -  (f^nf  (5) 

a 

where  j  is  the  current  density  maximum  and  j /a  is  the  voltage 
drop  (dU/dz  =  Ez)  or  the  electric  field,  one  obtains,  since 
Ez  is  constant  over  the  cross  sectional  area  of  the  throat  and 

A  A 

therefore  also  equals  to  Ez  =  j/cr, 

2ura  q  (  r  o )  =  nr02  j - -  I  Ex  (6) 

r  a  n  +  2 

A 

Herein  I  is  the  arc  current  and  o  is  the  maximum  value  of  the 
electrical  conductivity.  Since  for  fully  ionized  plasmas  a  can 
be  approximated  by 

a  a  v  T 
1  e 

where  Y  *  1.55  •  10"3  [ - 5 

V  m  KJ/ 

b  H.  0.  Schrade,  M.  Auweter-Kurtz  and  H.  L.  Kurtz,  "Stability 
Problems  in  magneto  Plasmadynamic  Arc  Thruster",  AIAA  18th  Fluid 
Dynamics  and  Plasmadynamics  and  Lasers  Conference,  AIAA-85-1633, 
Cincinnati,  OH,  July  1985 
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one  obtains  for  the  heat  flux  per  unit  length 


2*ro  qr(ro)  =  — - 


n  +  2  0  .65  •  103  I2 


and  if  one  replaces  the  maximum  electron  temperature  by 


T  =  257  (■—  +-A  ) 1//5  (— )2/5 

e  v  n  '  r  o 


it  is 


1  +  2 
n 


2nr0  qr(r0)  “(-^-4)3/10 


5.0 


i  o~ 2  (—  )7//5  r- 
ro' 


[- 

Lm 


(7) 


Q 
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The  results  of  tne  heat  load  measurements  on  the  innermost  seg¬ 
ment  of  the  nozzle  throat  (see  Fig.  14)  of  the  IRS  DT2  thruster 
witn  ro  =  1.2  cm  were  now  compared  with  that  obtained  by  the 
above  equation  (see  Fig.  14).  The  agreement  is  very  good  at 
higher  currents  (~  4000  A)  and  deviates  by  a  factor  of  about  two 
at  low  currents  (~  1000  A).  At  low  currents  the  arc  channel  may 
not  fill  the  entire  throat  such  that  the  radially  transported 
heat  from  the  arc  is  partially  carried  away  along  the  axis  by  a 
cold  gas  layer.  With  increasing  currents,  however,  the  arc  cross 
section  assumes  more  and  more  that  of  the  throat,  while  the  cold 
gas  layer  diminishes  and  the  radial  heat  losses  of  the  arc  col¬ 
umn  agree  with  the  measured  heat  load  to  the  inner  throat  walls. 


H.O.  Schrade,  M.  Auweter-Kurtz  and  H.L.  Kurtz,  "Basic 
Processes  of  Plasma  Propulsion",  Final  Scientific  Report,  AFOSR 
ftrant  82-0298,  Jan.  1987 


Fig.  14  Comparison  between  calculated  and  measured  radial  heat 
flux  per  unit  axis  length  within  the  nozzle  throat. 


26 


A  reevaluation  of  Hiigels  onset  explanation6  has  been  conducted 
by  means  of  the  flow  discharge  calculations  described  above.  The 
radial  particle  density  was  determined  by  equalizing  the  radial 
Lorentz  forces  with  the  pressure  gradient,  and  for  a  given  mass 
flow  rate  m  the  arc  current  was  increased  until  the  pressure  or 
particle  density  at  the  innermost  anode  point  dropped  to  zero. 
Plots  of  constant  particle  density  lines  for  which  this  "plasma 
starvation"  at  the  anode  occurs  are  shown  in  Fig.  15a  and  b  for 
m  =  0.0006  kg/s  with  a  critical  current  of  I  =  3600  A  and 
for  m  =  0.0015  kg/s  with  a  critical  current  of  6020  A,  respec¬ 
tively.  Comparison  between  several  of  these  calculations  with 
measured  onset  conditions'^  show  excellent  agreement  (see  Fig. 
16)  . 


6  H.Hiigel,  "Zur  Funktionsweise  der  Anode  im  Eigenfeldbeschleu- 
niger",  DFVLR  Report  FB-80-20,  1980 

'  H.  L.  Kurtz,  M.  Auweter-Kur tz  and  H.  0.  Schrade,  "Self  Field 
MPD  Thruster  Design  -  Experimental  and  Theoretical  Investiga¬ 
tions",  18th  international  Electric  Propulsion  Conference,  AIAA 
85-2002,  Alexandria,  VA,  1985 
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3 . 2  Cylindrical  Thruster 

In  addition  to  the  steady  state  nozzle  type  thruster  calcula¬ 
tions,  a  time  dependent,  fully  two-dimensional,  cylindrical  MPD 
channel  flow  calculation  has  been  conducted.  This  numerical 
treatment,  which  assumes  rotational  symmetry,  friction  free 
flow,  no  azimuthal  current  or  velocity  and  the  electric  field 
perpendicular  to  the  lengthwise  limited  coaxial  electrodes  (see 
Fig.  17)  has  been  described  in  Appendix  III.  Besides  the  example 
already  shown  there,  another  one  with  a  shorter  current  rise 
time  (~  10  us  instead  of  200  us)  and  a  lower  mass  flow  rate 
(m  =  5  g/s  instead  of  20  g/s)  has  been  conducted  and  is  presen¬ 
ted  here. 

Based  on  the  time  dependent  current  profile  as  shown  in  Fig.  18, 
the  current  contour  lines  and  the  pertinent  plasma  properties 
like  pressure,  temperature  and  velocity  were  calculated  as  a 
function  of  r  and  z  for  argon  as  propellant.  At  begin  (t  =  0) 
the  current  was  assumed  strictly  radial  and  homogeneously  dis¬ 
tributed  across  the  electrode  surfaces  and  the  axial  flow  field 
correspondingly  homogeneously  distributed  across  the  cylindrical 
tube  cross  section.  The  initial  temperature  and  the  Mach  number 
of  the  flow  were  taken  T(t=0)  =  10000°K  and  Ma(t=0)  =  1.5,  re¬ 
spectively.  In  Fig.  17  the  current  contour  lines  are  shown  after 
60  us;  the  temperature,  pressure  and  density  fields  after  10  us 
and  after  60  us  each  are  plotted  in  Figs.  19,  20  and  21,  respec¬ 
tively.  The  axial  velocity  vz,  the  Mach  number  and  the  mag¬ 
netic  induction  field  are  given  in  Figs.  22,  23  and  24,  respec¬ 
tively. 

Interesting  to  note  is  the  fact  that  the  current  contour  lines 
show  a  compression  not  only  on  the  downstream  side  but  faintly 
also  on  the  upstream  side  of  the  cathode.  This  effect  could  not 
be  observed  in  the  previous  example  (Appendix  III).  The  tempera¬ 
ture  and  the  pressure  field  show  peaks  on  both  ends  of  the  cath¬ 
ode  during  the  current  rise  period.  These  peaks,  however,  flat¬ 
ten  out  with  time  indicating  more  likely  a  certain  insufficiency 
in  the  initial  conditions  and/or  in  the  computer  code  than  in  a 


distribution  within 


netic  induction  distribution  within  the  channel  at 
10  ys,  b)  60  ys. 


i 


1 


'1 


-  37  - 


Fig.  25:  Thrust  versus  time. 
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Fig.  28  Localization  of  the  different  positions  for  the  time 
dependent  considerations  (Figs.  25,  26,  27). 


real  effect.  Nevertheless,  the  computations  so  far  yield  reason¬ 
able  overall  results  as  shown  by  the  thrust  versus  time  curves 
in  Fig.  25  and  by  the  temperature  and  pressure  versus  time  plots 
in  Figs.  26  and  27.  in  these  last  three  figures  the  differences 
in  temperature  and  in  pressure  at  the  middle  position  and  the 
near  anode  position  are  negligible;  therefore,  they  lie  for  both 
positions  on  one  curve,  respectively.  The  localization  of  the 
different  positions  is  shown  in  Fig.  28. 

The  thrust  as  plotted  in  Fig.  25  is  defined  by 

r 

Th  =  2n  j  (  pv  2  +  p)  r  dr 

r  z 
c 

It  is  interesting  to  note  that  after  about  40  ns  the  thrust  at 
the  end  of  the  channel  reaches  a  value  which  is  a  little  higher 
than  the  sum  of  the  gas  dynamic  and  electromagnetic  thrust. 
After  that  it  decreases  slowly  and  will  reach  this  value  (of  gas 
dynamic  and  electromagnetic  tnrust)  beyond  60  ns.  During  the 
current  rise  period  (0  -  10  ns)  the  thrust  at  the  channel  end 
does  not  change,  and  only  after  more  than  about  60  ns  does  the 
thrust  adjust  to  its  constant  steady  state  value. 

Similarly,  the  temperature  curve  (see  Fig.  26)  at  the  end  of  the 
channel,  which  is  practically  the  same  curve  for  all  three  posi¬ 
tions  7,8,9,  reaches  its  steady  state  constant  value  only  after 
more  than  60  ns.  The  temperatures  at  the  middle  positions  1,2,3 
of  the  electrodes  increase  at  first  to  about  16000°K  and  after 
that  drop  slowly  to  about  11000°  -  12000°K  after  again  more  than 
60  ns . 

The  corresponding  pressure  curves  in  Fig.  27  show  the  same  over¬ 
all  trend  as  the  temperature  curves.  The  pressure  at  the  channel 
end  reaches  29  mbar  after  about  55  ns. 

The  time  x  which  elapses  after  reaching  overall  steady  state 
conditions  con  now  be  estimated  by  the  length  of  the  active 
channel  portion  (given  by  210  mm  -  60  mm  =  0. 1 5  m)  divided  by 
the  average  speed  of  sound  (for  argon  at  about  12000°K, 
a  *  2»103  m/s)  as 
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t  =  =  7. 5  •  10“ b  s  =  75  (iS 

2-10 

This  time  is  in  good  agreement  with  the  adjustment  times  as  cal¬ 
culated  for  the  thrust,  temperature  and  pressure  according  to 
the  two-dimensional,  time  dependent  MPD  channel  flow  calcula¬ 
tions  presented. 
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4.  List  of  reports  and  papers  prepared  during  research  period 
(1  August  1986  -  31  July  1987) 

Tne  following  papers  and  student  theses  were  prepared  which  di¬ 
rectly  relate  to  the  research  work  of  Grant  APOSR  86-0337. 

Student  Theses  (Diplom) 

P.  Lehninger,  "Designing  and  Building  of  a  Test  Pig  in  Order 
to  Measure  Cathode  Erosion  Caused  by  an  Arc  Pulse".  IRS, 
Sept.  1986 

A.  Isselhorst,  "Calculation  of  the  Expansion  Process  in  a 
Plasma  Thruster  by  Accounting  for  Electromagnetic  Compres¬ 
sion".  IRS,  May  1987 

Papers 


M.  Auweter-Kurtz ,  H.L.  Kurtz,  H.O.  Schrade  and  P.C.  Sleziona, 
"Numerical  Modeling  of  the  Flow  Discharge  in  MPD  Thrusters", 
19th  International  Electric  Propulsion  Conference,  AIAA-87- 
1091,  May  1987 

H.  O.  Schrade,  M.  Auweter-Kurtz  and  H.  L,  Kurtz,  "Cathode 
Phenomnena  in  Plasma  Thrusters",  19th  International  Electric 
Propulsion  Conference,  AIAA-87-1 096 ,  May  1987 
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6. 1  Appendix  I 

Calculation  of 

where  e 


the  Dimensionless  Quantity  e 


8  n 


-  J  p  vz2  dAz  -  [l  +  f  ( n)  ] 


m  qI  A 


(1.0) 


r 
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Fig.  II  Illustration  of  the  current  density  lines  and  the 
Axial  Velocity  within  a  Spot  Discharge 


For  the  high  pressure  plasma  and  vapor  region  within  the  crater 
(see  dashed  control  volume  in  Fig.  1.1)  the  equation  of  motion 
follows  by 

7(pvv  +  p)  =  ]  x  B  (1.1) 


P2a.i-w*i&«x  _ 


T 


•w 
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and  after  integration  over  the  control  volume  V  and  applying 
Gauss'  theorem,  one  obtains 

/  {pvv  +  i>}  dA  =  f  J  x  B  dV  (1.2) 

A  V 

The  surface  integral  on  the  left  side  can  now  be  split  into  two 
terms,  the  integral  over  the  cross  sectional  area  of  the  dis¬ 
charge  channel  at  the  crater  orifice  (Az  =  itrQ2)  and  the  sur¬ 
face  integral  over  the  inner  spherical  crater  surface  (Ac;  see 
Fig.  1.1).  For  the  z  component  of  eq.  (1.2)  one  therefore 
obtains 


J  {pvv  +  p }  dA  u  +  J  {pvv  +P}dAu  =  J  j  x  B  dV*u  (1.3) 

A  Z  Z  A  2  2  V  2 

z  c 

where  u  is  the  unit  vector  in  the  z  direction.  If  one  now  de- 
z 

notes  the  vapor  and  plasma  pressure  of  the  inner  crater  surface 
by  pc  and  the  pressure  and  the  axial  velocity  at  the  orifice 
by  pz  and  vz,  respectively,  one  can  replace  eq.  (1.3)  by 


i 

i 

1 

1 

j 


"  /  {pv  2  +  p  )  dA  +  /  {p  }'(di  -u  )  =  /  3  x  $  dV*u  (1.4) 

AZZZACCZV  z 

z  c 

Herein  the  negation  sign  of  the  first  integral  stems  from  d^  u  = 

z  z 

-dAz;  one  neglects  friction  effects  and  assumes  that  at  the 
inner  crater  surface  the  average  randomized,  thermal  speed  of 
the  particles  prevails  over  their  convective  velocity.  Since  now 

dA„*u„  =  +  dA  =  +  r  drd9  (1.5) 

c  z  z 

where  9  is  the  azimuth  angle  of  the  cylinder  coordinates  r,9,z, 
it  follows  for  the  z  component  of  the  force  balance  eq.  (1.3) 

/  pv  2  dA  =  /  ( P  -  PJ  dA  -  /  3  x  B  dVu  (1.6) 

A  Z  ZA  C  2  2  V  Z 

z  z 

The  pressure  pc  is  now  taken  as  constant,  and  by  accounting 
for  the  magnetic  or  pinch  pressure  within  the  cross  sectional 
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area  Az  and  by  denoting  the  ambient  pressure  through  p®, 
one  obtains 


2 

/  pvz  dAz  =  (pc  -  pj. Az  - 


2  r  +  ■+■  ,  > 

I  -  /  3  x  B  dV'u, 


(1.7) 


V 


In  order  to  calculate  the  volume  integral  over  the  electromag¬ 
netic  forces  in  eq.  (1.7),  one  now  assumes  that  the  current  den¬ 
sity  is  evenly  distributed  over  the  inner  crater  surface.  The 
shape  of  the  crater  shall  be  approximated  by  a  spherical  segment 
which  can  be  characterized  by  the  angle  <J>  o  and  the  radius 
ro  =  rcsin<)>o  (see  Fig.  1.2).  By  introducing  the  Maxwell's 
stress  tensor  ^Maxwell/  one  now  obtains  for  the  volume 
integral 


/  3  x  B  dV 
V 


TMaxwell  ‘dA 


2  (J.  o 


/  { 2BB  -  B2  1} 


dA 


(1.8) 


where  the  magnetic  induction  vector  B  is  given  by  the  sum  of  the 
self  field  3  and  the  applied  field  So 

O 

B  =  §s  +  B0  (1.9) 


For  the  integrand,  one  may  write 

2BB  -  B2  T  =  2B^Bs  +  2EToBq  +  2B^B0  +  2B_oS 

-  (Bg  +  Bo2  +  2Bs*B0)  1 


s 


(1.10) 


Since  now  Bo  can  be  considered  constant  and  the  vector  Bs  is 
normal  to  dA  over  the  entire  surface  of  the  control  volume,  the 
integral  (1.8)  reduces  to 

/  3  x  B  dV  =  — -  /  { 2bTb  o  -  (B  2  +  2B  B0)  l}  *dA  (1.11) 

H  0  ®  o  o 


and  the  z  component  of  this  force  to 


H.  0.  Schrade,  "Magnetoplasmadynamic  Effects  in  Electric 
Arcs",  Interim  Scientific  Report,  Grant  AFOSR  82-0298,  August 
1983,  p.  1 09f f . 
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/  3  x  B  dV  u 


-351  /  V  “*•«%> 


(I. 12) 


The  other  terms  vanish  since  B  and  the  divergence  of  B  are 

s  z  s 

zero.  Splitting  the  surface  integral  in  two  parts  as  before,  one 
obtains 


/  J  X  5  av-u2  -  +  I  /  6. 2  aaz  -  /  «s2  <«*<.•".>  1 

z  c 


(1.13) 


where 


dAz  =  r  dr  d0  and 


(dA  *u  )  =  r  sin<j>  cos4>  d 4>  d0 
c  z  c 


(a) 
(1.14) 

(b) 


For  a  straight  channel  with  a  current  density  distribution  given 
by  a  paraboloid  of  grade  n,  one  obtains  for  the  integral  over 
the  cross  sectional  area  of  the  crater  orifice 


^ —  /  B  2  dA, 

2nu  .  s  3 

z 


-  TTZ  I  f(n) 


(1.15) 


with 


_L  { ( n  ♦  2 )  2  -  + 


(1.16) 


For  a  flat  current  density  profile  with  n  =  10,  f(10)  =  0.655; 
for  a  steep  profile  with  n  =  2,  f ( 2 )  -  0.917. 

The  second  surface  integral  in  eq.  (1.13)  follows  by 


1  r  ^  2u4>q 

27 7T  /  Bs2  (dA  u.)  =  /  /  Bs2  sin4>  cos*  d*  d0 

2n0  s  c  z  0=o  0=Q  s 

(1.17) 

Since  now  the  current  density  is  evenly  distributed  over  the  in¬ 
ner  spherical  crater  surface  it  follows  according  to  Maxwell's 


B  =  _ 
s  2n  r  sm<|> 
c 


(1.16) 


where 


2 7C  $0 

!{$)  =  /  /  j  r  sin$  d 4>  d9  =  2n  r  2  j 

0  0  c  c 

=  2 nr c 2  j  [l  -  cos  $  ] 

Since  the  total  current  carried  by  the  spot  is 
I  =  2nrc2  j  [  1  -  cos<t>o] 


o 

/  sin$  d <t> 

0 


(1.19) 


(1.20) 


one  obtains 

I  (  $  )  =  I  | 


cos<{> 
cost  0 


(1.21) 


and  according  to  (1.18) 

no  I  1  1  -  cos4> 

3  =  _  _  _  _ 

s  2 it  rc  sin<t>  1  -  cos<t>o 


(1.22) 


Introduced  into  eq.  (1.17)  yields  according  to  Standard  Integral 
Tables  9 


1 

2no 


<d;Vz) 


Ho  2  2 

011  ( 1  —cos  <t>  0  )  2 


/ 

<P  =  0 


*  0  (  1-COS(>>)  2 

sini)> 


cosi)  d<j> 


HO  2  4 

011  (1-COS<t>o)2 


4>o 

(1+2  ln(cos  y~ ) 


cos<|>  o 


HO 


2CUo) 


(1.23) 


The  dimensionless  function  C <  o )  is  plotted  in  Fig.  1.3  for  the 
range  0  <  $o  <  n/2.  The  limit  angle  4>o  is  about  zero  for  a  flat 
crater  and  for  a  hemispherical  crater  it/2,  c ( 4> o )  varies  therfore 
between  0.250  and  0.114  depending  on  the  shape  of  the  crater. 


e.g.  Samuel  M.  Selby,  "Standard  Mathematical  Tables",  16th 
Edition,  The  Chemical  Rubber  Co.,  Cleveland,  OH,  1968 
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The  volume  integral  over  the  electromagnetic  forces  in  the  z-di- 
rection  follows  now  according  to  eqs .  (1.13),  (1.15)  and  (1.23) 
by 

H  0 

/  3  x  B  dV-£  =  q-t  I  2  {f  {  n)  -  2  C  ( 4>  o )  }  (1.24) 

V  z  Bu 

This  introduced  into  eq.  (1.7)  yields 

0 

/  pvz2  dAz  =  (pc  -  Poj)A2  -  8^  I2  {1  +  f(n)  -  2  C  ( ♦  o )  }  (1.25) 

A 

z 

and  hence  the  dimensionless  quantity  e  [eq.  (1.0)]  can  be  writ¬ 
ten 


Pc  ~  P' 

HO 

si  1  i 

pc  -  p 
HU 

81  1  J 


2  [1  +  f(n)  -  C(<Do)  ] 


3.3 


(1.26) 


(1.27) 


The  last  approximation  is  well  justified,  since  f(n)  =  0.8 

(±  0.2)  and  C(<t>o)  =  0.15  (±  0.04).  The  factor  3.3  may  therefore 
vary  by  about  ±  15%  depending  on  n,  the  grade  of  the  current 
density  paraboloid,  and  on  $o,  which  characterizes  the  crater 
shape . 
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Abstract 

Processes  at  the  arc  cathode  attachment  decisively 
determine  the  entire  discharge  behavior  of  almost 
all  arc  devices  and  therefore  also  of  MPD  and/or 
arc  jet  thrusters.  Cne  well  known  process  occur ing 
on  spotty  arc  attachments  in  a  transverse  magnetic 
field  is  the  fact  that  the  cathode  spots  move  or 
jimp  in  tne  direction  opposite  to  the  Lorentzian 
rule.  In  pulsed  thruster  devices  with  cold  cathodes 
and  very  likely  also  in  continuously  running 
thrusters  with  so-called  thermionic,  seemingly  dif¬ 
fuse  attachments  of  hot  surfaces,  the  arc  attach¬ 
ment  consists  of  many  high  current  density  spots. 
These  spots  can  stick  or  spread  upstream  and  there¬ 
by  overheat  the  insulating  material  of  the  back- 
plate  of  the  thruster.  In  this  paper  an  explanation 
of  the  phenomena  of  spot  motion  is  presented. 

Introduction 

Within  pulsed  gas-fed  plasma  thrusters  with  coaxial 
electrodes,  one  often  observes  that  relatively  high 
erosion  takes  place  at  the  upstream  side  there  the 
center  cathode  penetrates  the  backplate  insulator 
(see  Fig.  1). 

This  fact,  together  with  the  arc  traces  left  on  the 
cathode  surface,  indicate  that  the  arc  cathode  at¬ 
tachment  which  consists  of  many  highly  instationary 
hot  spots  sticks  preferably  to  the  cathode  base. 
Scanning  electron  micrographs  of  cathode  surfaces 
which  were  charged  by  high  current  pulses  of  up  to 
10  kA  reveal  typical  craters  of  about  1  im  to  al¬ 
most  100  urn  in  diameter,  which  confirms  the  exist¬ 
ence  of  many  seemingly  statistically  distributed 
hot  spots.  Any  one  of  these  active  spots  may  be 
characterized  by  a  crater-like  pit  with  an  inner 
molten  layer,  as  shown  in  Fig.  2. 

Within  the  pit  exists  a  high  pressure  plaana  caused 
by  ohmic  heating  and  heavy  evaporation  of  the  elec¬ 
tron  emitting  inner  surface  layer.  This  high 


Fig.  1  Illustration  of  the  area  with  relatively 
high  erosion. 
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pressure  plaana  expands  out  into  the  ambient  lower 
pressure  gas  sphere  and  forms  a  plasma  jet.  Depend¬ 
ing  on  the  lifetime  of  the  spot  and  the  size  of  the 
molten  layer,  droplets  can  also  be  ejected  from 
such  an  active  site,  drastically  increasing  the 
surface  damage  and  the  erosion  rate  of  such  a  spot¬ 
ty  arc  attachment.  Since  the  ejected  plasma  is 
electrically  conducting,  one  must  also  consider 
that  a  current-carrying  channel  oaanates  from  the 
crater  orifice.  The  current  density  of  the  orifice 
can  assume  values  of  tp  to  the  order  of  10  h/m 
and  possibly  even  higher1,  and  the  power  density 
within  an  active  spot  correspondingly  reaches  val¬ 
ues  of  to  about  10 16  kw/nr  within  a  region  of  a 
few  micrometers. 

The  arc  attachment  consists  now  of  many  of  these 
spots  jmping  seemingly  statistically  over  the 
cathode  surface.  As  shown  experimentally  in  many 
different  arc  devices,  the  first  time  to  Stark  and 
since  then  by  nuserous  researchers,  ,  ,  the  spotty 
arc  attachment  moves  in  the  direction  opposite  to 
Ampere's  law  if  an  applied  magnetic  field  is  pres¬ 
ent  and  is  oriented  parallel  to  the  cathode  sur¬ 
face.  Since  its  discovery,  this  well  known  but 
still  not  well  understood  phenomenon  of  arc  retro¬ 
grade  motion  has  sparked  the  imagination  of  many 
scientists,  stirring  them  to  develop  various  hy¬ 
potheses.  -iI  A  presentation  and  scrutiny  of  all 
these  interpretations  would  extend  beyond  the  scope 
of  this  paper,  and  a  good  many  have  been  experimen¬ 
tally  refuted  and/or  do  not  rest  on  solid  physical 
foundations.  Even  the  many  experimental  findings  at 
first  did  not  fully  correspond  with  each  other; 
however,  since  the  scanning  electron  microscope  for 
the  examination  of  the  arc  footprints  on  the  cath¬ 
ode  surface  is  available,  the  experimental  results 
become  unequivocal  but  reveal  even  more  complexity 
as  originally  assured. 

An  additional  phenomenon  was  observed  first  by 
Smith12  and  later  by  Kesaev13  in  a  diverging 
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Fig.  2  Scheme  of  a  spot  discharge  with  a  straight 
discharge  channel. 


magnetic  field  and  afterwards  was  experimentally 
proven  to  occur  al9o  in  an  inclined  constant  mag¬ 
netic  field.  ,  It  shows  that  if  the  magnetic 
field  lines  are  somewhat  inclined  with  respect  to 
the  cathode  surface,  the  spots  discontinuously 
drift  away  from  the  retrograde  direction.  This 
drift  motion  was  extensively  investigated  by  Robson 
and  is  therefore  named  after  him. 

In  this  paper  an  analysis  is  presented  which 
uniquely  explains  these  strange  phenomena.  The  con¬ 
sequences  an  the  spot  behavior  in  a  coaxial  plaaaa 
thruster  are  discussed. 

Qualitative  Explanation 

The  vapor  and  plasma  jet  which  proceeds  from  a 
superheated  cathode  spot  is  electrically  conductive 
and  thus  contains  the  intrinsic,  current-carrying 
channel;  the  latter  electrically  connects  the  cath¬ 
ode  spot  and  thus  the  cathode  with  the  plasna  of 
the  interelectrode  space  and  thereby  with  the 
anode. 

First  let  us  consider  an  axially  symmetrical  dis¬ 
charge  channel  which  converges  to  the  electrode, 
having  its  origin  in  the  cathode  spot.  The  magnetic 
field  lines  of  such  a  current-carrying  channel  are 
the  concentric  circles  aroind  the  channel  axis  and, 
in  accordance  with  the  cowerging  lines  of  the  cur¬ 
rent  density  vector  fiel induce  an  electromag¬ 
netic  volume  force  ([]  x  BJ  force),  which  has  com¬ 
ponents  normal  and/or  parallel  to  the  axis.  If  one 
plots  the  values  of  these  forces  acting  radially 
inward  on  a  cross  sectional  plane,  then  one  obtains 
a  volcano-like  profile,  as  seen  in  Fig.  3. 

Oh  the  basis  of  this  force  profile,  one  obtains  a 
pressure  increase  in  the  plasma  channel  towards  the 
discharge  center,  which  is  known  as  “pinch  pres¬ 
sure”.  The  average  pressure  increase  in  the  channel 
is  calculated  as 
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^Hag  ‘  Ij  (1) 

where  J  represents  the  current  density  averaged 
over  the  cross-sectional  area  and  I  is  the  electric 
current  flowing  through  the  channel  (40  »  4x*10~ 
Vs/Am  magnetic  permeability).  Actually,  this  is  a 
well  known  effect;  it  shows,  however,  that  in  low 
pressure  discharges  the  magnetic  pressure  increase 
can  be  of  a  much  greater  magnitude  than  the  sur¬ 
rounding  pressure,  and  that  self-magnetic  or  rather 
plasmadynamic  effects  dominate  over  electrostatic 
effects  in  the  quasi-neutral  discharge  channel. 
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!  swmetric,  straight  spot  dis- 
x  B  force  distribution. 


Fig^  4  ^Curved  asymmetric  spot  discharge  with  radi¬ 
al  J  x  B  force  distribution. 

If  this  channel  is  curved  (see  Fig.  4),  then  the 
magnetic  field  lines  on  the  concave  side  of  the 
channel  are  condensed  and  on  the  convex  side  are 
spread  out.  This  means,  however,  that  the  radial 
[]  x  Bj  forces  on  the  concave  side  become  larger 
than  on  the  convex  side,  and  the  original,  axially 
symetrical,  volcano-like  force  profile,  distrib¬ 
uted  over  the  cross-sectional  area  of  the  channel, 
becomes  asymnetrioal.  This  asymmetry  has  the  effect 
that  the  cathode  vapor  plasma,  under  high  pressure 
and  ohmically  heated,  can  now  no  longer  be  held  to¬ 
gether  by  the  magnetic  forces;  it  expands  or 
streams  out  on  the  weak-foroe  side  of  the  channel 
and  is  thereby  diverted  off  the  axial  direction. 
Analogous  to  the  operation  of  rocket  propulsion,  in 
which  the  repulsion  of  the  outflowing  gw  produces 
surface  forces  an  the  walls  of  the  oortustion  cham¬ 
ber  and  propells  it  into  the  opposite  direction  of 
the  flow,  the  discharge  channel  is  now  driven  in 
the  direction  opposite  to  the  deflected  vapor  and 
plasma  beam,  and  the  chmaber  walls  of  the  rocket 
^ngige  may  be  compared  with  the  confining  radial 
]  1  B  force  configuration,  l.e.  the  axial  vapor- 
arid  plasma  jet  flowing  out  from  the  cathode  spot  is 
deflected  in  the  direction  of  mailer  force  fields, 
while  the  intrinsic  current-carrying  channel  is  de¬ 
flected  in  the  opposite  direction  (see  Pig.  4). 

under  certain  conditions,  this  deflection  of  the 
current-carrying  channel  can  become  so  strong,  that 
it  contacts  and  heats  up  the  cathode  surface  at  a 
new  point  neighboring  the  spot.  Cathode  material  is 
thereby  vaporized  also  at  this  point,  and  ultimate¬ 
ly  a  new  electrode  contact  in  the  form  of  a  new 
cathode  spot  arises  (see  Fig.  5). 

In  fact,  an  electrically  conducting  plaama-v*»r 
mixture  exists  outside  of  the  intrinsic  current- 
carrying  channel,  so  that  a  large  part,  but  not  all 
of  the  current  carried  by  the  spot  will  flow 
through  the  intrinsic  current  channel.  The  current 
density  in  the  intrinsic  charnel  is  always  much 
larger,  however,  than  that  outside  of  the  channel 
considered. 


Fig.  3  Rotational 
charge  with  radial 


Fig.  5  Spot  propagation  due  to  bending  of  the 
current  channel. 

Fixed  Coordinate  System 

Fig.  7  Illustration  of  the  channel  bound  (x,y,z) 
and  the  fixed  (X£,»yL'2L)  coordinate  systems 
with  respect  to  the  laboratory  frame  of  reference. 
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Fig.  6  Straight  spot  in  a  transverse  magnetic 
field.  Plaana  jet  is  deflected  in  the  amperian  di¬ 
rection. 

If  the  discharge  channel  is  situated  in  an  external 
magnetic  field  running  parallel  to  the  cathode  sur¬ 
face,  the  radial  []  x  5]  force  profile  can  become 
asymmetrical  without  additional  curving  of  the  dis¬ 
charge  channel  (see  Fig.  6).  The  effect  of  this 
asymmetry,  however,  as  in  the  case  discussed  above, 
is  a  deflection  of  the  vapor  and  plasma  bean  in  the 
amperian  direction,  while  a  resulting  gas  dynamic 
surface  force  acts  on  the  current-carrying  channel 
in  the  opposite  (retrograde)  direction.  It  can  now 
be  shown  that  under  certain  conditions  this  gas  dy- 
namic  force  in  the  retrograde  direction  predomi¬ 
nates  over  the  lorentz  force  (I  x  B0)  in  the  emper- 
ian  direction.  The  channel  axis  is  thereby  curved 
in  the  retrograde  direction  until  the  discharge 
channel  ultimately,  as  in  the  previous  exanple 
(shown  in  Fig.  5),  contacts  and  heats  up  the  cath¬ 
ode  surface.  In  this  way,  a  new,  distinct  cathode 
attachment  point  arises  on  the  retrograde  side  of 
the  original  spot.  The  cathode  attachment  point  re¬ 
produces  itself  accordingly  in  the  retrograde  di¬ 
rection. 

Theory 

One  considers  any  arbitrary  configuration  of  a 
current-carrying  channel  within  a  fixed  laboratory 
systa*  (xL.yL«*L:  •«  Pig-  7)-  The  theory 
starts  frcm  the  equation  of  motion,  the  continuity 
equation  and  Harwell's  equation 


3  »  electric  current  density 

B  *  magnetic  induction 

Che  now  replaces  the  velocity  vector  v,  of  the 
fixed  system  (xL,yL,zL)  by 

vL  *  v  +  vs  (6) 

i.e.  by  the  sum  of  the  plasma  velocity  v,  measured 
within  the  channel  bougd  system  (x,y,z;  see  Fig.  7) 
and  the  axis  velocity  v  which  has  only  transverse 
components  with  respectsto  the  channel  axis 
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After  integration  of  the  equation  of  motion  ever  a 
channel  segment  of  length  As  (see  Figs.  7  and  8) 
one  eventually  obtains  the  force  balance  foe  the 
channel  segment  in  the  form 

dv  + 

(*>  -  i  *  (8) 

where  tn  •  /  p  dv  is  the  mass  of  the 

AV  plants  within  the 

senwnt  of  volume  AV 
and  length  As 
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Pig.  10  Plot  of  the  functions  f(n)  and  g(n). 


Fig.  8  Curved  discharge  channel  segment  of  length 
4s  =  R4*  (R  *  ralius  of  curvature  vector). 


X  4s  is  the  sun  of  all  forces  which  act  on  the 
segment  and  accounts  for 

a)  all  gas  dynamic  surface  forces  like  propelling, 
friction,  and  pressure  effects 

b)  the  volume  forces  due  to  electromagnetic  effects 
between  the  current-carrying  elements  inside  and 
outside  the  segment  and 

c)  all  fictitious  forces  like  Ooriolis  and  centrif¬ 
ugal  forces  due  to  an  axial  flow  within  a  curved 
and/or  twisted  current-carrying  plasma  channel. 

This  normal  force  ii  per  unit  channel  length  has 
now  been  calculated1®  by  modeling  the  current  den¬ 
sity  distribution  within  the  channel  by  any  "eccen¬ 
tric  paraboloid  of  grade  n*.  Such  a  distribution 
follows  from  a  concentric  paraboloid  through  shift¬ 
ing  the  maximum  to  a  point  of  the  channel  cross 
section  with  the  coordinates  ( r»,  8m)  as  shown 

in  Fig.  9. 

This  modeling  allows  a  very  good  approximation  of 
any  real  current  density  distribution  by  adjusting 
the  grade  n  and  the  eccentricity  coordinates  (rn, 
0M) .  This  is  the  more  justified  since  in  the 


final  expression  of  the  force  Xi  the  current  den¬ 
sity  j  shows  cp  only  within  the  integral  expres¬ 
sions 

ih  J,  %2  -  is 12  f(n>  f1  +  r 

(9) 
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Fig.  9  Current  distribution  modeled  by  a)  a  con¬ 
centric  and  b)  an  eccentric  paraboloid  of  grade  n. 


■i  i  r  r  r  icose,  _  i  ‘m 

I  *  I  .  !  -1  r0  sin8^  r  r  ~  n+3  ro  sine 


where  is  the  self  magnetic  field  induced  by  the 
diannel  current  1  and  which  depends  on  the  current 
density  distribution  within  the  channel;  rdrdS  « 
dAj  is  a  surface  element  of  any  cross  sectional 
area  of  the  channel;  a,p  are  measures  of  the  eccen¬ 
tricity  of  the  current  density  (moments  of  j  in  the 
x~  and  y-  direction);  and  f(n)  and  g{n)  are  func¬ 
tions  of  the  grade  n  ^nd  are  plotted  in  Pig.  10. 

The  transverse  force  Xj^  follows  finally  in  the  form 


^  1*0  1 0 

X1  ’  1  I1  +  f(n)  +JT  (f(n)9(n)(n*3)  +  0® 

+  7  +  k  W  p  ar  (?*®0)  ■  J  ^z 


+  R  3s  J  p  vz(vy  ♦  V  dA*>  3 
\  K 

*  (fear  f  ** 

♦  Ras  ^ p  v*  <vx ♦  vsx>  <**>-!  -  jt  ■ 

(»i) 


*z  *  wo  channel  cross  section 
dA^  dro 

31“  “  2*r°  ar  channel  divergence  or  convergence 
2*  ro 

/  ...  dA  ■  /  /  ...  r  dr  d8  sumaation  over 


8"0  r>0 


the  cross-sec¬ 
tional  area 


S’ro>ghi  A«>t 


gg  *  channel  twist,  change  of  the  twist  angle  along 
the  channel  axis 

v  ,v  ,v  »  velocity  components  measured  within  the 
‘  channel  bound  coordinate  system 

v  ,  v  »  transverse  velocity  components  of  the 
”  channel  axis 


I  *  I  u  *  current  vector  in  z -direction 
(along  axis) 

$u  »  applied  magnetic  field  vector 

This  force  £j_  per  unit  channel  length  is  now  among 
others  a  function  of  the  local  channel  geometry 
like  channel  radius  m,  axis  curvature  R,  channel 
twist,  convergence  or  divergence  and  also  a  func¬ 
tion  of  the  channel  orientation  with  respect  to  an 
applied  magnetic  field  Bo  and/or  flow  field. 


In  this  calculation  one  assumes  that,  inside  the 
discharge  channel  and  in  its  immediate  vicinity, 
the  time  derivation  of  the  channel  bound  impulse 
density  can  be  neglected,  or  that  flow  and  pressure 
£iel<l  adjust  in  a  relaxation-free  manner  to  the 
]iB  forces  and  vice  versa,  this  assumption,  how¬ 
ever,  is  valid  as  long  as  the  change  of  the  channel 
dianeter  j(2ro)/dt  or  the  motion  of  the  discharge 
channel  |v  I  is  small  compared  to  the  average  speed 
of  sound  within  the  considered  channel  segment, 
i.e. 


d(2r o) 
dt 


«  a 


(12) 


Hie  validity  of  this  inequality  can  be  proven  from 
case  to  case  and,  as  discussed  elsevhere, 1 6  is  ful¬ 
filled  for  even  extreme  conditions. 

Based  on  the  knowledge  of  X i  one  can  now  determine 
■balanced"  and  "stable"  discharge  channel  configur¬ 
ations. 

Fbr  a  balanced  current-carrying  glasma  channel  con¬ 
figuration  the  transverse  force  Xj_  along  the  dis¬ 
charge  axis  must  be  zero,  i.e.  the  gas  dynamic  sur¬ 
face  forces  and  all  body  forces  like  electromagnet¬ 
ic  field  and  fictitious  forces  must  be  mutually 
balanced.  Such  a  channel  configuration  for  instance 
has,  according  to  eq.  (11),  a  straight  axis  with 
R  »  Considering  now  the  discharge  channel  that 
onanates  from  a  fixed  spot  crater  (see  Fig.  11), 
one  can  define  its  orientation  with  respect  to  an 
applied  magnetic  field  Bo  by  the  inclination  angle 
5.  If  the  channel  axis  becomes  slightly  disturbed 
in  any  direction  *  by  a  small  deviation  £  off  its 
straight  configuration,  the  axis  becomes  curved  and 
the  deviation  or  disturbance  vector  £  points  now 
always  in  the  direction  of  the  radius  of  curvature 
vector  R.  Disturbances  by  lateral  motion  of  the 
electron  emitting  crater  with  respect  to  the  dis¬ 
charge  channel  are  disregarded,  because  any  such 
motion  is  physically  unrealistic,  since  the  dis¬ 
charge  channel  must  be  fixed  to  the  electron  emit¬ 
ting  hot  spot.  Nevertheless,  many  explanations  of 
arc  motion  are  based  on  such  a  dubious  lateral 
effect. 


By  definition,  such  a  balanced  discharge  channel 
will  be  called  "stable"  if  any^disturbance  £  initi¬ 
ates  a  transverse  force  X  »  X(£)  *  0  which  acts 
against  6  and  drives  the  channel  back  in  its  origi¬ 
nal  balanced  configuration.  The  mathematical  re¬ 
quirements  for  a  balanced  and  stable  discharge 
channel  configuration  follow  therefore  by 

4  •  t  <  0  (13) 


Fig.  11  Illustration  of  a  disturbed  current-carry¬ 
ing  channel  at  the  cathode  surface  with  inclined 
magnetic  field. 

If  this  requirement  is  not  met,  $  acts  in  the  di¬ 
rection  of  £  and  forces  the  channel  to  deviate  more 
and  more  from  its  balanced  state. 

Let  us  now  discuss  by  means  of  this  stability  theo¬ 
ry  the  dynamics  of  a  spot  discharge. 

Spot  Motion 

With  the  introduction  of  the  inclination  angle  l 
and  the  disturbance  apgle  *  as  shown  in  Fig.  11, 
the  transverse  force  jq  which  acts  on  the  current- 
carrying  plasma  channel  of  a  spot  discharge  can  be 
replaced  by 
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(15) 


k  »  k(n)  «  (n+3)f(n)g(n)  +  1 


(16) 
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_I _ ,  avg.  ylf  magnetic  Induction  field 

roBo  applied  magnetic  induction  field 


Herein  one  assumes  that  any  disturbance  t  leads 
only  to  a  bent  (see  Fig.  11)  but  not  to  a  twisted 
current-carrying  channel.  According  to  eqs.  (11) 
and  (14),  a  straight  channel  (with  R  *  «)  is  again 
a  balanced  one,  since  Xi  •  0.  As  mentioned  before, 
the  direction  of  any  disturbance  £  points  now  in 
the  R  direction;  therefore,  a  balanced,  straight 
channel  is  also  stable  if  the  braced  expression  of 
the  R/R  component  of  Xi  is  larger  than  aero.  Neg¬ 
lecting  the  higher  order  terns  of  ro/R,  one  there¬ 
fore  obtains  as  the  stability  requirement 

«  -  £  cosC  cos*  >  0  (17) 

Cne  may  now  distinguish  between  the  following  three 
cases 
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1.  No  applied  magnetic  field,  |B0j  ~  ^  »  0 


stable  by  bending  more  and  more  in  the  retrograde 
direction. 


2.  So  parallel  to  the  cathode  surface,  5  *  0  . 

3.  So  inclined  with  respect  to  the  cathode  surface, 
5*0. 

Case  1;  J  B0 1  ~  ^  «  0 

In  this  case  the  stability  requirement  reduces  to 

t  >  0  (stable)  (18) 

which  is  identical  with  that  of  the  MPD  arc  in  sec¬ 
tion  3.2.  Here  the  stability  criterion  does  not  de¬ 
pend  on  the  direction  of  the  disturbance  given  by  * 
(see  Pig.  11).  Within  eq.  (15)  one  may  now  relate 
the  average  impulse  transport 


along  the  channel  axis  with  the  y vapor  and  plaaaa 
pressure  py  of  the  spot  crater1  and  the  aibient 
pressure  p».  Therefore,  the  stability  requirement 
(18)  can  also  be  expressed  in  the  form 

MO 

Pv  -  P-  >  3-3  Ti  la  3s  (19) 

where  I  is  the  current  carried  by  the  spot  and  5S 
is  the  Iwerage  current  density  at  the  crater  ori¬ 
fice;  the  factor  3.3  is  valid  for  a  hemispherical 
crater  and  changes  to  about  3.0  for  a  flat  crater. 

The  pressure  difference  between  the  vapor  pressure 
in  the  spot  and  the  anbient  pressure  must  be  there¬ 
fore  more  than  about  three  times  larger  than  the 
average  magnetic  pressure  (a,  *  uu/8*  I  j  )  in  or¬ 
der  to  maintain  a  straight,  stable  spot  discharge. 
If  this  requirement  is  no  longer  fulfilled  or  if 
c  <  0,  the  channel  becomes  unstable.  Any  snail  dis¬ 
turbance  will  cause  the  channel  to  bend  more  and 
more,  and  eventually  it  will  contact  the  cathode 
surface  at  a  nearby  site  off  the  original  spot  (see 
Fig.  5).  The  new  site  is  heated  cp  and  a  new  spot 
is  created.  Again,  for  this  new  spot,  the  stability 
requirement  (18)  or  (19)  must  be  fulfilled  in  order 
to  maintain  a  balanced,  stable  spot.  Since  in  this 
case  the  stability  criterion  does  not  depend  on  the 
direction  of  the  disturbance,  and  singe  one  may  as¬ 
sure  that  all  possible  directions  of  6  are  equally 
probable,  the  spots  should  move  or  jump  statisti¬ 
cally  over  the  cathode  surface. 

Case  2;  5q  parallel  to  the  cathode  surface,  5*0. 


Por  a  <  0,  the  current  density  max  inns  is  shifted 
away  from  the  center  of  curvature  maximim  towards 
the  convex  side  of  the  channel.  According  to  cri¬ 
teria  (20)  and  eq.  (21),  the  channel  should  now 
bend  preferably  in  the  osperian  direction.  Che  can, 
however,  show  that  by  accounting  for  the  higher 
order  term  ( r0/R) ak  in  eq.  (14),  the  channel  this 
time  assures  a  curved,  stable  configuration.  The 
possible  consequence  of  this  behavior  shows  that 
the  channel  does  not  come  in  contact  with  the  cath¬ 
ode  surface  on  the  amperian  side  of  the  spot.  This 
is  not  the  case  for  a  retrograde  disturbance  with 
a  >  0,  for  which  the  channel  is  unstable  and  bends 
more  and  more.  Cbservable  current  channels  should 
therefore  be  amperically  bent  or  straight  while  the 
observation  of  a  retrograde  channel  configuration 
would  require  high  time  resolutions. 

In  conclusion  to  Case  2,  a  spot  prefers  to  jugp  or 
move  either  in  the  retrograde  or  the  amperian  di¬ 
rection,  depending  on  vhether  the  current  density 
maximun  is  shifted  towards  the  concave  or  the  con¬ 
vex  side  of  the  discharge  channel.  Because  of  the 
unstable  channel  behavior  during  retrograde  motion 
and  the  stable,  bent  channel  behavior  during  amper¬ 
ian  motion,  one  may  conclude  that  the  average  jimp 
distance  is  analler  and  the  nrnber  of  jmps  more 
frequent  for  the  retrograde  moving  spots  than  for 
the  amperian  moving  ones.  Pb^  large  magnetic  fields 
(XJc  «  1,  which  means  about  Bo  »  1  Tesla  for  an  ex¬ 
perimental  value1  of  1/ro  *  10 '  Vm)  the  mechan¬ 
isms  for  oqperian  and  retrograde  spot  motion  are 
practically  identical  and  therefore  in  both  direc¬ 
tions  equally  probable. 

Case  3:  So  inclined  with  respect  to  the  cathode 
surface,  5*0. 

According  to  eq.  (20)  the  current  channel  is  pref¬ 
erably  unstable  with  respect  to  a  disturbance  in 
the  retrograde  direction  and  assures  a  bent  but 
stable  configuration  for  an  amperian  disturbance. 
The  principle  difference^  non  is  given  through  the 
force  component  in  the  lit  direction,  which  in 
the  first  order  of  Cq/R  follows  by 


In  this  case,  the  stability  criterion  follows  by 

e  -  “  cos*  >  0  (20) 

i.e.  the  channel  stability  depends  now  on  the  di¬ 
rection  of  the  disturbance  (*).  moreover,  the  force 
component  in  the  I  x  R  direction  of  eq.  (14)  be¬ 
comes  unequal  to  aero  in  the  first  order.  It  is 

M  t2  ,  f_jt 

X1  (in  fx£)  *  "  U*  IT  V  8ln*  TO  (21) 


and  which  turns  the  disturbed  or  curved  channel 
around  its  original  straight  axis  either  in  the 
retrograde  (a  >  0)  or  amperian  direction  (a  <  0). 
tow,  however,  the  disturbance  6  finally  does  not 
point  exactly  in  these  directions  but  deviates  or 
drifts  awey  fro*  that.  The  drift  angle  follows  from 
the  requirement  that  for  this  position  the  turning 
force  in  the  t  x  R  direction  must  be  aero.  Prom 
eq.  (22)  follows  therefore 

i  *ro 

drift  ‘  -  IS  ST  ^  <23> 


tor  a  >  0,  i.e.  the  current  density  max lmm  is 
shitted  towards  the  center  of  curvature,  the  chan¬ 
nel  is  less  likely  to  be  stable  for  any  disturbance 
towards  the  retrograde  aide  -n/2  <  *  <  tx/2  than 
for  a  disturbance  towards  the  aeperian  side  t/2  < 
♦  <  3»/2  (or  -x/2  >  ♦  >  -3x/2)  which  can  be  veri¬ 
fied  by  means  of  criteria*  (20).  for  any  disturbed 
or  curved  channel,  there  acts  now  a  force  in  the  ♦ 
♦ 

I  s  R  direction  [eq.  (14)]  such  that  the  curved 
charnel  always  turns  into  the  retrograde  direction, 
i.e.  i  points  finally  in  the  opposite  I  x  Bo  direc¬ 
tion.  The  channel  therefore  becomes  preferably  un¬ 


for  a  spot  action  towards  the  retrograde  side  it  is 
e  >  0,  and  since  the  channel  in  the  direction  of 
the  current  is  supposed  to  converge  (drg/ds  <  0), 
the  drift  angle  4drtft  is  proportional  to  5  as 
long  as  the  angles  are  mull.  The  retrograde  moving 
spots  drift,  therefore,  inder  a  —11  angle  with 
respect  to  the  proper  retrograde  direction  towards 
the  right  (in  the  direction  of  Be)  if  the  Bo  field 
points  upward,  and  drifts  towards  the  left  (oppo¬ 
site  to  Be)  if  the  field  vector  points  downward 
(a—  Pig.  12a,b).  This  spot  behavior  has  been  ax- 
periman tally  observes  and  as  mentioned  before  is 
known  —  fobaon  drift  motion.  , 15 
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Fig.  12  Retrograde  spot  motion  with  an  inclined  So 
field. 

Discussion 

Let  us  consider  the  cathode  attachment  of  an  MPD 
arc,  consisting  of  many  active  spots  which  are 
somehow  distributed  over  the  cathode  surface,  as 
shown  for  instance  in  Fig.  13.  Any  single  spot  is 
now  affected  by  all  the  other  ones  such  that  a  mag¬ 
netic  field  acts  on  it.  Assume  now  that  all  the 
spots  are  statistically  homogeneously  distributed 
aroind  and  along  the  cathode  surface.  In  this  case 
the  magnetic  induction  field  is  azimuthal  and  de¬ 
pends  on  how  much  current  is  drawn  by  the  spots 
downstream  of  the  considered  spot. 


Cathode  Rode 


Fig.  13  Scheme  of  the  spotty  arc  cathode  attach¬ 
ment. 

For  a  cathode  rod  of  length  1  and  diameter  d*. 
covered  statistically  by  many  spots,  the  current 
through  the  rod  at  any  distance  z  from  the  back- 
plate  follows  by 


where  1  is  the  total  current,  d^  is  the  cathode 
diameter,  1  ia  the  length  of  the  rod  and  z  is  the 
distance  from  the  backplate. 

therefore  the  azimuth  ail  magnetic  induction  field 
for  a  spot  at  a  distance  z  fra*  the  backplate  is 
about 


Ihis  has  the  consequence  that  Bo  decreases  down¬ 
stream.  For  any  spot,  therefore,  one  has  to  take 
into  account  a  transverse  magnetic  field,  the  con¬ 
sequence  of  such  a  field,  as  shown  in  the  previous 
section  under  ’Case  2’,  is  a  preferred  notion  of 
the  spots  in  the  retrograde  or  upstrean  direction. 
With  increasing  magnetic  field  the  quotient  k  in 
eq.  (17)  decreases,  which  leads  to  an  enhancement 
of  the  instability  behavior  and  therefore  to  a 
faster  spot  motion.  Since  the  crater  size  of  fast 
moving  spots  can  be  expected  to  be  smaller  than 
that  of  slowly  moving  ones,  the  spot  size  should 
increase  towards  the  cathode  tip.  This  has  been  ob¬ 
served1  in  seemingly  diffuse  cathode  attachments 
and  indicates  that  possibly  all  attachment  mechan¬ 
isms  are  spotty  and  locally  highly  unstable. 70 

In  sunmarizing  the  consequences  on  plasma  thruster 
cathodes  the  following  should  be  noted. 

>  A  cathode  spot  is  the  more  stable  and  sticks 
preferably  to  that  area  where  the  vapor  pressure 
of  the  cathode  material  is  high.  Since  insulator 
material  outgases  or  evaporates  easier  than  the 
thoriated  tungsten  rod  cathode,  the  arc  attach¬ 
ment  usually  sticks  to  the  boundary  cathode  in¬ 
sulator,  or  cathode  backplate. 

*  Within  a  magnetic  field  parallel  to  the  cathode 
surface,  an  individual  spot  jumps  preferably  in 
the  retrograde  direction,  i.e.  upstream. 

•  The  jump  distances  of  the  retrograde  moving  spots 
are  smaller  than  that  of  the  amperian  moving 
ones. 


•  If  the  magnetic  field  lines  are  inclined  with  re¬ 
spect  to  the  cathode  surface,  the  spots  drift 
away  frcm  the  retrograde  direction. 

•  The  upstream  boundary  of  a  seemingly  diffuse  arc 
attachment  depends  on  the  conditions  under  which 
retrograde  arc  spot  motion  ceases.  Those  condi¬ 
tions  are  not  fully  known  today. 

•  While  the  current-carrying  channel  of  a  spot  dis¬ 
charge  in  a  transverse  magnetic  field  may  bend 
and  eventually  move  d isoont inuously  in  the  retro¬ 
grade  direction,  the  plaaaa  jet  is  deflected  in 
the  anperian  direction  and  may  cause  a  bright 
layer  along  the  cathode  rod.  In  this  layer,  evap¬ 
orated  cathode  material  should  be  present. 
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Abstract 

Two  theoretical  models  for  calculating  the  current 
and  flow  distributions  in  self-field  MPD  thrusters 
are  beinj  developed  and  will  be  applied  to  evaluate 
the  effects  of  geometry,  propellant  type,  scale  and 
other  parameters  on  the  thruster  performance. 

1)  For  aontinusus  thrusters,  a  stationary  oode  is 
being  developed.  The  extended  Ohm's  law  is  used 
to  calculate  the  current  contour  lines,  are!  a 
one-dimensional,  two-ccmponent  expansion  flow 
model  is  used  to  obtain  the  velocity,  tempera¬ 
ture  and  pressure  distributions  for  calculating 
the  gas  properties,  trfuch  are  again  used  in 
Ohm's  law.  The  differential  equation  is  solved 
by  means  of  a  finite  difference  method  for  tne 
geometry  of  the  nozzle-type  plasma  thruster  OT2, 
whicn  has  been  investigated  at  IRS  in  a  steady 
state  as  well  as  in  a  quasi-steady  state  mode. 

An  integration  over  the  volume  and  thermal 
focces  equals  to  the  tnrust.  'hie  calculated  cur¬ 
rent  density  distribution  and  the  computed 
thrust  are  compared  with  experimental  results. 

2)  for  tne  starting  phase  of  the  steady  state  MPD 
thrusters  as  well  as  for  pulsed  thrusters,  a 
time  dependent,  fully  two-dimensional  code  is 
being  developed.  It  uses  a  modified  McCormack  FD 
method  in  cylindrical  coordinates  to  calculate 
the  time  dependent  flow,  temperature  and  pres¬ 
sure  fields.  The  current  and  magnetic  contour 
lines,  vhich  must  be  known  for  the  energy  and 
impulse  input  in  this  code,  are  calculated  with 
an  instationary  Ohm’s  law  solved  in  a  manner 
similar  to  the  computer  oode  for  continuous 
thrusters.  This  oode  is  used  for  a  simplified 
cylindrical  thruster  model. 

for  both  methods  the  results  are  discussed. 

1 .  Introduction 

The  increase  in  onboard  electrical  power  in  the 
near  future  is  again  raising  interest  in  self-field 
mpd  thrusters.  Despite  their  simplicity  in  design 
and  power  conditioning,  they  are  handicapped  even 
today  by  the  shortcomings  of  low  efficiency  and  a 
performance  limit  known  as  "onset"  which  restricts 
the  flow  rate  at  a  given  current  value.  Therefore, 
since  the  beginning  of  the  design  and  investigation 
of  self-field  MPD  thrusters  at  the  end  of  the  six¬ 
ties,  effort  has  been  made  to  unterstand  the  phys¬ 
ics  of  the  MPD  discharge  and  the  development  of  ap¬ 
propriate  codes.  But  the  MPD  problems  cannot  be 
solved  easily,  because  in  contrast  to  the  hardware 
simplicity  of  these  devices,  the  physics  involved 
are  of  estrone  complexity.  Tb  mention  only  a  few: 
real  gas  effects,  I.e.  cold  propellant  is  heated  up 
in  an  arc  and  the  gas  has  to  be  dissociated  and 
ionised ;  the  problems  related  with  the  electrodes 
are  only  partly  knowi  or  solved;  the  plans  is  In  a 
thermal  non-equilibrium;  because  of  the  low  densi¬ 
ties  and  high  temperatures  the  Reynold's  numbers 
are  low  m)  the  friction  cannot  be  neglected;  the 
magnetic  Reynold's  number  Is  too  high  to  separate 
flow  and  discharge;  etc.  this  list  is  by  no  means  a 
ample te  one.  Added  to  these  physical  complexities 
are  numerical  difficulties:  the  necessary  rolling 
of  different  partial  differential  equation  systoas. 


rumely  elliptical  and  hyperbolic  ones;  non-lineari¬ 
ties;  steep  gradients;  and  ao  on.  And  let  us  not 
forget  difficult  geometries.  Therefore,  due  to  the 
great  complexity  of  the  physical  conditions,  for 
all  codes  established  to  solve  at  least  parts  of 
these  problems,  many  and  great  simplifications  must 
be  made.  The  resulting  codes  include  for  example 
one-dimensional1,',  quasi-two-^ imensional  and  sim¬ 
plified  two-dimensional  ones1-  . 

In  this  paper,  two  further  approaches  are  present¬ 
ed:  first,  a  partly  two-dimensional  one.  The  sim¬ 
plified,  two-dimensional  electromagnetic  calcula¬ 
tion  follows  front  Ohm's  law,  similar  to  Ref.  5,  and 
the  flow  is  modeled  by  a  one-dimensional,  non-equi¬ 
librium  expansion.  The  equations  are  solved  only 
for  the  steady  phase.  The  second  approach  uses  com¬ 
parable  electromagnetic  equations,  but  models  the 
flow  actually  two-dimensionally  with  a  modified 
MacCormack  oode.  Both  approaches  neglect  electrode 
effects  and  friction. 


Steady-State  MPD  Thruster  Investigations 


In  order  to  predict  the  overall  performance  of  con¬ 
tinuously  ruining,  self-field  MPD  thrusters,  a 
semi-two-dimensional  model  calculation  and  computer 
progran  has  been  developed.8.  In  a  first  step  this 
program  has  been  used  to  investigate  the  current 
contour  lines  and  the  thrust  within  the  nozzle-type 
thruster  OT2  (Fig.  2.1),  uhich  is  being  developed 
at  this  institute. 


2. 1  Current  Contour  Lines 

In  this  approach,  the  current  contour  lines  are 
calculated  in  r-  and  z-  dependence  while  the  flow 
field  is  taken  as  a  one-dimensional,  frictionless 
expansion  flow,  assuming  isothermal  behavior  for 
the  electron  component  and  lsentropic  expansion  for 
the  heavy  particles  like  ions  and  neutrals.  This 
assumption  is  justified,  since  primarily  the  elec¬ 
tron  component  is  heated  by  Chmic  energy  input,  and 
since  within  the  nozzle  expansion  flow  the  energy 
exchange  between  the  lightweight  electrons  and  the 


Fig.  2.)  Self-field  nozzle-type  MFD  thruster  DT2 
(capacity  about  5  kA  or  15  N  thrust  with  argon). 
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2.10 


heavy  particles  like  ions  and  neutrals  is  fairly 
weak.  Therefore,  reaction  collisions  (ionization 
and  recombination  processes,  etc.)  are  neglected. 
Hence,  the  model  approach  distinguishes  between 
electron  and  heavy  particle  temperatures  and  as¬ 
sumes  a  frozen  flow  from  the  nozzle  throat  on  down¬ 
stream. 

Within  the  nozzle  throat,  one  assures  a  "fully  de¬ 
veloped*  arc  flow  regime,  i.e.  the  gas  and  plaana 
velocity  and  the  current  density  vector  have  only 
an  axial  component  and  the  arc  oolunn  is  defined  by 
only  radial  heat  transfer  (and  possibly  by  radial 
radiation)  losses.  Moreover,  from  the  nozzle  throat 
on  upstream,  one  no  longer  distinguishes  between 
electron  and  heavy  particle  temperatures,  since 
within  this  higher  pressure  regime  a  tight  coupling 
between  all  components  is  more  likely  to  be  estab¬ 
lished  than  in  the  low  pressure  downstream  nozzle 
regime.  In  spite  of  the  neglect  of  the  radial  de¬ 
pendency  of  the  flow  and  species  tanperatures,  how¬ 
ever,  this  model  approach  turns  out  to  be  an  effec¬ 
tive  tool  for  predicting  the  current  contour  lines 
within  and  the  thrust  and  specific  impulse  of  a 
nozzle-type  self-field  MPD  thruster. 

Assumptions: 

•  steady  state  conditions 

•  rotational  symmetry 

•  no  azimuthal  current 

•  singly  ionized,  quasi-neutral  plaana 

•  electric  current  density  normal  to  electrode 
surfaces 

•  from  cathode  to  nozzle  throat:  thermal  equilib¬ 
rist 

•  within  the  nozzle:  frozen  flow,  ambipolar  ex¬ 
pansion  flow  without  friction,  electrons:  iso¬ 
thermal,  ions:  adiabatic 

•  outside  tne  nozzle:  hyperbolic  expansion  flow, 
electrons  and  ions  adiabatic 
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where  f  is  a  slowly  varying  function  depending  on 
current  distribution. 

The  last  three  equations  are  derived  in  detail 
under  the  above  assumptions  in  itef.  8. 

Rewriting  Ohm's  Law  by  means  of  Maxwell's  equa¬ 
tions,  one  obtains  a  veqtor  equation  for  the  mag¬ 
netic  induction  vector  B  in  the  form 
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For  a  stream  function 
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with  respect  to  the  rotational  symmetry  and  the 
zero  azimuthal  current,  the  elliptic,  partial  dif¬ 
ferential  equation  of  2nd  order  follows  from  eq. 
2.11  in  the  form 
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Basic  Equations: 
Maxwell 's  eqs. : 
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Ohm's  law:  .  _  . 

i  -  -  a  +  v  x  si  -  f  [j  *  i] 

2.5 

Equation  of  State: 

p-  pfi.  ti  ♦  (fr’i 

2.6 

The  function  ?(r,z)  *  const,  represents  now  a  cur¬ 
rent  contour  line,  since  B  =  B0  ~  I(r)/r  and  I(r) 
is  the  electric  current  carried  through  a  cross 
sectional  area  of  mr.  The  proper  boundary  condi¬ 
tions  for  T  follow  frco  the  geometry  of  the  thrust¬ 
er  walls  and  electrodes. 

A  solution  of  the  equation  system  2. 1  -  2.13  re¬ 
quires  as  input  data: 

•  Geometry  of  Thruster  A(z) 

•  Mass  Plow  Rate  m 

•  Electric  current  I 

and  yields  the  following  results: 

•  Current  Contour  Lines 

•  Thrust ,  $ 
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Specific  impulse,  Isp 


Expansion  Relation: 


Step  1:  By  setting  the  integral 
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Nozzle  Throat  Conditions 
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in  the  expansion  relation  eq.  2.8  to  aero, 
one  obtains  together  with  the  nozzle 
throat  conditions  2.9  and  2.10  a  first  ap¬ 
proximation  of  v(s),  p(z),  p(t)  and 
nj(z)  -  ne(z)  within  the  nozzle.  Out¬ 
side  the  nozzle  one  assuaes  a  hyperbolic 
expansion  flow,  which  am  ana  vr  *  0 
(Fig.  2.2). 

With  the  knowledge  of  the  flow  and  plamaa 
properties  (in  first  approximation) ,  o,  9 
and  v  are  alao  knoxxi,  as  well  am  their 
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Pig.  2.2  Scheme  of  the  density  model. 

Inside  the  nozzle:  one-dimensional. 

Outside  the  nozzle:  elliptical. 

derivatives.  With  these  quantities  one  now 
solves  the  differential  equation  2.13  by 
means  of  a  finite  differential  method 
(Gauss-Seidel)  for  the  current  contour 
lines  in  first  approximation.  In  Fig.  2.3 
the  current  contour  lines  are  shown  for 
different  currents  between  1.3  and  4  kA. 
According  to  Maxwell's  equation,  one  also 
c^btayis  frcm  ?  the  current  density  and  the 
]  x  B  force  (configuration  within  the  dis¬ 
charge. 

Step  2:  Che  solves  the  expansion  jelajion  (eq. 

2.8)  nunerically  by  taking  (]  x  B)  within 
within  the  integral  frcm  Step  1.  With  that 
the  velocity  v(z)  and  the  quantities  p(z), 
p(z),  ripiz)  are  calculated  in  second  or¬ 
der  approximation  and  from  that  the  new 
quantities  for  a,  8  and  v  with  their  de¬ 
rivatives. 

With  these  new  quantities,  one  now  solves 
the  differential  equation  2.13  for  the 
current  contour  lines  by  means  of  the  sane 
miner! cal  method  as  in  Step  1.  Repeating 
Step  2  yields  an  iterative  solution  of  the 
equation  system  2.1  -  2.10. 

The  calculation  shows  that  after  two  steps 
the  current  contour  lines  with  I  -  4  kA  do 
not  change  anymore  and  even  a  first  step 
approach  gives  good  results.  FOr  lower 
current  rates  the  deviation  is  even 
analler. 

A  comparison  between  calculated  and  measured  cur¬ 
rent  pattern  (Fig.  2.4)  shows  a  fairly  good  agree¬ 
ment.  The  fact  that  in  the  measured  case  the  anode 
attachment  is  more  concentrated  downstream  of  the 
nozzle  end  than  in  the  calculated  case  may  be  at¬ 
tributed  to  the  theoretical  model  which  neglects 
the  boundary  layer  along  the  nozzle  and  anode 
walls. 

2.2  Thrust 

According  to  the  known  g  x  S  force  configuration 
and  the  flow  conditions,  one  can  now  calculate  the 
thrust. 

The  thrust  of  a  self-field  MFD  accelerator  is  com¬ 
posed  of  two  effects;  first  the  bub  of  all  gas- 
dynanic  surface  forces  on  the  inner  plenum  chamber 
and  on  the  inner  nozzle  wall,  and  second  the  sun  of 
all  voluae  effects  given  by  the  magnetic  field 
forces  which  act  between  the  current-carrying  plas¬ 
ma  and  the  entire  power  circuit,  including  the  co¬ 
axial  electrodes.  Hence  it  is 
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Fig.  2.3  Current  contour  lines  for  different  cur¬ 
rents  with  a  mass  flow  of  0.8  g/s. 

where  {pvv  +  p)  is  the  gas  dynamic  force  per  wit 
area;  Ag  represents  the  surface  of  all  internal 
walls  like  plenim  chamber  with  the  gas  Inlet  cross 
sections,  the  cathode  surface  and  the  inner  solid 
wall  of  the  nozzle;  V  is  the  current-carrying tvol^ 
use.  The  second  term  is  the  sub  over  all  j  x  B 
forces  which  act  between  the  current-carrying  solid 
leads  of  the  electric  power  circuit  and  the  cur¬ 
rent-carrying  plaaaa;  in  the  second  equation  po  is 
the  chancer  pressure,  A^  the  nozzle  throat  area 
and  ojh  the  thrust  coefficient  Matching  the  gas 
dynamic  effects.  It  is  by  definition 
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Fig.  2.4  a)  Calculated  and  b)  measured  current 
contour  lines  of  the  MPD  discharge  in  the  nozzle 
type  OT2-IRS  thruster  for  I  =  1,3  kA  and  m  -  0.8 
g/s  argon. 

CT  =  7rx~  J  (tvv  +  p}*dA  2.15 
1  Po\  A 

S 

The  calculated  thrust,  based  on  the  current  density 
and  magnetic  field  determined  with  the  procedure  as 
described  in  section  2.1,  is  now  compared  with  the 
experimental  results  in  Fig.  2.5. 


Fig.  2.5  CDnpariaon  between  calculated  and  meas¬ 
ured  thrust  of  the  OT2-IRS  thruster.  T^^ax  is 
the  maximum  possible  electromagnetic  thrust. 

They  show  excellent  agreement  foe  op  ■  1  with  the 
experimental  measurements  in  the  upper  current 
range  (I  >  3000  A),  and  even  down  to  1300  A  the  de¬ 
viation  is  less  than  30%.  The  pure  electromagnetic 
thrust  is  given  by  the  lower  curve  in  Fig.  2.5.  Che 
can  see  that  the  electromagnetic  thrust  reaches 
about  50%  of  the  total  thrust  at  a  current  of  about 
4500  A. 


3.  Time  Dependent,  Two-Dimensional  >ffp  Thruster 
Investigations 

In  addition  to  the  steady  state  MFD  theory,  a  real 
two-d^ensional  numerical  code  has  been  devel¬ 
oped.  This  code  al lows  the  calculation  of  the 
time  dependent  plaana  properties  during  the  start¬ 
ing  phase  of  steady  state  MFD  thrusters,  as  well  as 
for  pulsed  mode  thrusters.  The  numerical  solution 
is  done  in  a  manner  similar  to  that  of  Jacobi12  for 
laser  applications.  In  a  first  approach,  this  MPD 
code  was  applied  to  the  cylindrical  MFD  thruster 
ZT1  being  developed  at  this  institute.  Fig.  3.1 
shows  the  drawing  of  ZT1,  and  Fig.  3.2  shows  the 
corresponding  simplified  scheme  for  the  calcula¬ 
tion. 


Fig.  3.1  Cylindrical  continuous  MPD  thruster 
ZT1-IRS. 
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Bild  3.2  Simplified  schone  used  for  the  calcula¬ 
tion. 

The  geometrical  configuration  on  which  the  amputa¬ 
tion  is  based  was  chosen  as  close  as  possible  to 
the  experimental  device.  The  radii  of  cathode  and 
anode  are  taken  according  to  the  dimensions  of  the 
real  device;  however,  the  length  of  the  cathode  is 
assumed  to  be  of  the  same  length  as  the  anode  and 
it  matches  the  length  of  the  three  anode  se^sents 
(see  Fig.  3.1).  The  model  is  assumed  to  be  an 
infinite  tube  of  an  insulator  with  a  central 
insulating  rod,  only  the  electrodes  will  be  of  zero 
resistance. 

for  the  time  impendent  calculation  of  the  MFD  flow, 
two  independent  codes  were  correlated  to  determine 
flow  field,  current  and  magnetic  field  dis¬ 
tribution.  These  two  codes  are  connected  in  the 
following  manner:  for  a  given  flow  field,  the  cur¬ 
rent,  and  hence  the  magnetic  field  distribution, 
was  determined.  With  these  results  the  flow  field 
equations  were  integrated.  In  the  next  time  step, 
this  new  flow  field  and  the  new  distribution  of  the 
thermal  properties  were  taken  to  calculate  the  new 
electromagnetic  field  distribution,  and  so  an. 

As  input  data  the  two  codes  require: 

•  geometry  of  the  thruster 

•  mass  flow  rate 

•  electric  current  as  a  function  of  time 
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and  at  the  starting  point: 

•  temperature  field 

•  Mach  number. 

This  yields  the  following  tune  dependent  results: 

•  flow 

•  temperature 

•  pressure 

•  tech  number 

•  magnetic  field 

•  current  density  distribution  within  the  channel 

•  electrothermal  thrust. 

3.1  Electromagnetic  Code 

The  electromagnet ic  properties  are  obtained  in  a 
similar  manner  as  frat  the  steady  state  thruster 
code,  using  the  following  assumptions: 

•  rotational  syrmetry 

•  no  azimuthal  current 

•  electric  field  normal  to  the  electrode  surfaces. 
Che  solves  Ohm's  law  with  respect  to  Maxwell's 
equations;  this  means  one  must  solve  the  elliptic 
partial  differential  equation  2.13. 

The  boundary  conditions  are  as  follows: 

Oh  the  insulator  surface  and  at  the  boundary  of  the 
integration  area,  the  stream  function  V  equals  zero 
for  the  case  of  the  downstream  insulator  and  equals 
¥max  f°r  the  upstream  insulator.  At  the  elec¬ 
trodes  the  boundary  condition  follows  from  the  re¬ 
quirement  that  the  electric  field  E  must  be  perpen¬ 
dicular  to  the  surface.  This  means  for  a  cylindri¬ 
cal  channel  as  considered  here  that  at  the  elec¬ 
trode  surtaces  Ez  equals  zero.  Therefore,  the 
boundary  condition  follows  to 

IdT  T  B*  ft?  „  3  , 

- vr - -  —  =  0  3. 1 

cuor  it  r  nor  8z 

With  this  closed  ooundary  condition,  the  computa¬ 
tion  has  been  done  with  a  Gauss -Seidel  code  and 
yields  as  a  result  the  magnetic  field  and  the  cur¬ 
rent  density  distribution  within  the  channel. 

3.2  Flow  Field  Code 

Fbr  the  description  of  the  time-dependent, 
two-dimensional,  supersonic,  ohmically  heated  flow 
the  conservation  equations  for  mass,  impulse  and 
energy  were  used  with  the  following  assumptions ; 

•  friction  free  flow 

•  rotational  symmetry 

•  no  azimutnal  current  and  velocity 

•  supersonic  flow  within  the  channel 

The  conservation  theorem  formulation  was  chosen  in 
order  to  satisfy  the  Rankin-Hugoniot  shock  condi¬ 
tions  implicity;  hence  the  description  of  shock 
depth  and  shock  velocity  are  very  good.  These  con¬ 
siderations  yield  the  following  non-linear  hyperbo¬ 
lic  differential  equation  system  with  cylindrical 
coordinates 
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where  the  indices  I L,  If  ,  |j  indicate  the  par¬ 
tial  differentiation  with  respect  to  time  and  to 
the  r-  and  z-directions. 

Boundary  Conditions: 

upetreaa:  at  that  boundary  all  flow  parameters 

are  assisted  to  be  constant  and  must 
be  set. 


downstream:  here  the  flow  parameters  are  extrapo¬ 

lated  because  of  the  calculated  field 
parameters . 

at  the  side:  no  normal  velocity;  normal  differen¬ 
tiations  equal  zero. 

Initial  conditions: 

At  t  *  0  the  flow  parameters  equal 
the  undisturbed  parallel  flow; 
uniform  current  density; 
tech  number  and  temperature  must  be 
set. 

For  the  numerical  solution  a  modified  tecOozmack 
code  was  used;  this  is  am  explicit  code  of  second 
order  which  uses  the  shock  capturing  method.  It  al¬ 
lows  a  fairly  high  energy  input  into  the  flow  due 
to  olmic  heating  and  a  steep  rise  in  the  impulse 
due  to  ]  x  B  forces.  Because  of  instabilities  at 
the  starting  phase  caused  by  initial  and  boundary 
conditions  on  the  upstream  boundary,  the  initial 
conditions  will  be  restricted  for  the  moment  to 
values  which  result  in  cowling  numbers  analler  than 
1.5  in  the  middle  of  the  electrodes. 

for  the  example  presented,  the  calculation  is  based 
on  the  following  initial  values:  T(0)  *  8000  K, 
Ma(0)  =  1.2,  with  a  mass  flow  of  20  g/s  argon.  The 
current  as  a  function  of  time  is  taken  as  a  sine 
square  fmction  sin2[xt/2to]  rising  fran  l(t=0)  =  0 
to  I(to*200us)  »  12  kA  and  is  taken  as  constant 
from  there  on. 


With  this  input  data  the  computation  yields  the 
time  dependent  current  density  distribution,  as  il¬ 
lustrated  in  Fig.  3.4.  Starting  with  radially 
parallel  and  homogeneous  current  contour  lines  at 
t  «  0,  the  current  is  driven  more  and  more  down¬ 
stream — stronger  near  the  anode  than  at  the  cath¬ 
ode-ami  bulges  at  the  electrode  ends.  The  dotted 
lines  are  those  at  20  us  after  starting,  and  the 
dashed  and  solid  lines  are  for  160  us  and  200  us, 
respectively. 


Fig.  3.4  Calculated  current  contour  lines  in  a  cy¬ 
lindrical  MFD  thruater.  Dotted  lines  at  t  •  20  is, 
dashed  after  160  us  anS  solid  after  200  us. 
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Fig.  3.5  illustrates  the  temperature  distribution 
100  us  and  200  us  after  starting.  The  strong  rise 
in  temperature  near  the  electrode  is  caused  by 
ohmic  heating  due  to  the  relatively  high  current 
density  there.  The  disturbance  at  the  cathode  end 
due  to  discontinuous  boundary  conditions  decreases 
with  time.  This  fairly  marked  disturbance  at  start¬ 
ing  is  one  reason  for  choosing  a  high  initial  tem¬ 
perature  . 


Fig.  3.5  Temperature  distribution  within  the  chan¬ 
nel  at  a)  100  us,  b)  200  us. 

A  radial  shock  front  formed  at  the  end  of  the  elec¬ 
trodes  moves  upstream  while  lowering  its  intensity; 
this  is  demonstrated  in  Fig.  3.6. 


Fig.  3.6  Pressure  distribution  within  the  channel 
at  a)  100  us,  b)  200  us. 


The  axial  and  radial  velocity  distribution  at  200 
us  are  plotted  in  Fig.  3.7.  At  both  ends  of  the 
cathode  there  is  a  steep  rise  in  the  radial  veloc¬ 
ity  vr,  while  the  axial  velocity  reaches  its  max¬ 
imum  at  the  beginning  of  the  cathode.  The  axial 
velocity  decreases  towards  the  cathode  end,  with  a 
low  rise  frcm  this  point  downstream  as  well  as  to¬ 
wards  the  anode. 


Fig.  3.7  Velocity  distribution  at  200  us.  a)  Axial 
velocity  vz,  b)  radial  velocity  vr. 

The  relationship  of  the  Mach  number  distribution  to 
the  temperature  and  velocity  distribution  at  200  us 
is  shown  in  Fig.  3.8.  The  Mach  number  set  equal  to 
1.2  at  the  upstream  boundary  reaches  an  initial 
maximum  at  the  beginning  of  the  cathode,  rises 
again  while  reaching  an  interim  maximum  at  about 
1/3  of  the  cathode  length  and  rises  once  more  be¬ 
hind  the  electrodes.  The  minimum  at  the  end  of  the 
cathode  is  caused  by  the  temperature  instability  at 
that  point.  It  is  remarkable  that  the  Mach  number 
is  below  1  within  a  wide  range,  especially  before 
the  end  of  the  electrodes. 


Fig.  3.8  Mach  number  distribution  at  200  us. 

As  described  in  section  2.2,  the  thrust  is  composed 
of  an  MPD  part  and  a  gw  dynamic  part  (see  eq. 
2.14a).  The  pure  MPD  thrust  can  be  calculated  by  an 
integration  of  the  j  x  B  forces  within  the  entire 
current-carrying  volume;  another  possibility  is  to 
use  the  well-known  formula13 


where  r;/rc  is  the  anode  to  cathode  ratio.  Both 
have  been  done  for  the  example  described  here  with 
a  coincidence  of  within  1%. 

In  Pig.  3.9  the  pure  MFD  thrust  (dashed  line)  and 
the  total  thrust  (solid  line)  are  plotted  depending 
on  time.  The  total  calculated  thrust  at  400  us  lies 
within  a  3»  deviation  fran  the  expected  thrust. 
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Fig.  3.9  Thrust  versus  time.  Solid  line:  total 
thrust.  Dashed  line:  electromagnetic  thrust. 

Even  more  troublesome  is  the  fact  that  the  conser¬ 
vation  of  mass  is  only  valid  within  a  10%  range. 
Fig.  3.10  shows  the  mass  flow  versus  time  at  four 
different  channel  cross  sections. 
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Fig.  3. 10  Mass  flow  versus  time  at  different  chan¬ 
nel  cross  sections. 

These  deviations  from  the  mass  conservation  law  are 
an  inherent  attribute  of  the  MacCormack  code. 


Concluding  Remarks 


In  order  to  gain  a  more  profound  understanding  of 
the  fundamental  process  occur ing  in  MFD  thrusters 
and  in  order  to  eventually  end  up  with  reliable  de¬ 
sign  criteria  and  predict  thruster  performances  un¬ 
der  various  conditions,  the  two  codes  described 
within  this  paper  have  been  developed. 


Che  was  developed  for  the  steady  state  condition 
and  the  other  for  the  time  dependent  description  of 
MPD  thrusters.  The  first  oode  was  tested  with  the 
OT2-IRS  thruster  design,  and  it  produced  fairly 
good  results;  especially  the  current  density  dis¬ 
tribution  agrees  quite  well  with  the  measurements . 
The  time  dependent  code  tested  with  the  cylindrical 
thruster  design  of  the  ZT1-XRS  thruster  produced 
fairly  good  results  as  well.  Both  codes  ham  weak 
points,  for  exosple  the  boundary  conditions,  real 
qas  effects,  and  geometry  modeling,  especially  the 
second  code  was  weak  in  terms  of  mass  conservation. 


In  a  later  research  period  the  physical  models  must 
be  refined,  and  other  miner  ical  flow  field  calcula¬ 
tion  codes  will  be  tested  In  terms  of  their  quali¬ 
fications  for  mag ne topi asn adynamic  thrusters. 
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